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We report fragmentation cross sections measured at 0° for beams of 14 N, ie O, 20 Ne, and 24 Mg 
ions, at energies ranging from 290 MeV/nucleon to 1000 MeV/nucleon. Beams were incident on 
targets of C, CH2, Al, Cu, Sn, and Pb, with the C and CH2 target data used to obtain hydrogen- 
target cross sections. Using methods established in earlier work, cross sections obtained with both 
large- acceptance and small- acceptance detectors are extracted from the data and when necessary 
corrected for acceptance effects. The large-acceptance data yield cross sections for fragments with 
charges approximately half of the beam charge and above, with minimal corrections. Cross sections 
for lighter fragments are obtained from small-acceptance spectra, with more significant, model- 
dependent corrections that account for the fragment angular distributions. Results for both charge- 
changing and fragment production cross sections are compared to the predictions of the Los Alamos 
version of the Quark Gluon String Model (LAQGSM) as well as the NUCFRG2 and PHITS models. 
For all beams and targets, cross sections for fragments as light as He are compared to the models. 
Estimates of multiplicity-weighted helium production cross sections are obtained from the data and 
compared to PHITS and LAQGSM predictions. Summary statistics show that the level of agreement 
between data and predictions is slightly better for PHITS than for either NUCFRG2 or LAQGSM. 

PACS numbers: 25.75.-q, 25.70.Mn, 25.60.Dz, 24.10.Lx, 98.70.Sa 



I. INTRODUCTION 

The Galactic Cosmic Rays (GCR) contain a small per- 
centage of heavy ions that contribute substantially to the 
dose and dose equivalent received in spaceflight par- 
ticularly in deep space outside the geomagnetosphere. 
As NASA's future plans are likely to include extended 
human missions in deep space, these exposures take on 
increased importance. A detailed understanding of the 
transport of these ions through matter is needed, as crew 
will typically be inside a modestly-shielded spacecraft, in 
habitats (conceivably with relatively thick shielding), or, 
in the case of a Mars mission, shielded by a combination 
of the CO2 atmosphere and a habitat. Fragmentation 
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cross sections play a key role in transport calculations 
and the resulting estimates of dose and dose equivalent 
behind shielding. These estimates can be highly uncer- 
tain Q due to propagation of cross section uncertain- 
ties. This can, in turn, be a central factor in limiting 
mission duration or rendering certain mission scenarios 
unfeasible. It is therefore essential that an accurate and 
precise database of the nuclear interaction cross sections 
is available to modelers for both code development and 
validation purposes. 

Historically, the space radiation community has fo- 
cused on heavy ions such as iron. This is understand- 
able, since the heavier GCR ions - iron in particular - 
contribute substantially to the dose and dose equivalent 
in unshielded deep space. However, from the standpoint 
of model completeness and reliability, it is also important 
that the fragmentation of lighter ions be well understood, 
since B, C, N, O, Ne, and Mg ions are much more abun- 
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dant than the heavier ions. Furthermore, beams of C 
d, IH and Ne ions @ have also been used in cancer ther- 
apy, where fragmentation plays a key role in limiting the 
(desired) localization of the dose to the tumor volume. 
We present cross sections for beams of 14 N ions at 290 
and 400 MeV/nucleon; 16 at 290, 400, 600, and 1000 
MeV/nucleon; 20 Ne at 290, 400, and 600 MeV/nucleon; 
and 24 Mg at 400 MeV/nucleon. Charge-changing (some- 
times referred to as "total" charge-changing) cross sec- 
tions and fragment production cross sections have been 
extracted from energy-loss spectra measured with sili- 
con detectors. The fragment production cross sections 
are partial cross sections. In 0° experiments such as 
this, each event has a "leading" (highest-charged) frag- 
ment. Absent angular acceptance corrections, the sum of 
these leading-fragment cross sections equals the charge- 
changing cross section for any given combination of beam 
ion and target. 

As in ourpreviously-reported fragmentation cross sec- 
tion data |6l— TlOlj . the charge-changing cross sections in 
the 250 to 1000 MeV/nucleon range are reproduced with 
reasonable accuracy by geometric models that treat the 
nuclei as overlapping spheres. In some cases, slight en- 
ergy dependence is observed in the data, and most of 
the models attempt to account for this, but on the whole 
it is a small effect and energy-independent models may 
suffice. In contrast, model predictions of the fragment 
production cross sections are, in general, not nearly as ac- 
curate as those for charge-changing cross sections. Here 
too there are subtle dependences on beam energy, and 
these are typically not well reproduced by the models. 
Older modelsje.g., Nilsen et al. [HI, NUCFRG2 [13], 
and EPAX2 [13| ) approximate fragment cross sections as 
monotonically decreasing functions of the charge change 
AZ (= Zbeam — -^frag), where Z is the nuclear charge. 
In many instances, this is fairly accurate in an average 
sense, over a limited range of AZ, but this approach 
cannot hope to reproduce important details seen in the 
data. Features missing in these older models include the 
enhanced production of even-Z fragment species relative 
to odd-Z species, suppression of F (Z = 9) production, 
and increases in cross sections for large AZ's. (For an ex- 
ample of these discrepancies, see Figure 15 of Ref. ficj.) 

Other high-quality data in the literature (e.g., fl4T4l6l |) 
report only fragment cross sections for (approximately) 
Zbeam/2. In contrast, in this and other articles we re- 
port the small charge change cross sections and also take 
the extra step to extract light-fragment production cross 
sections. This is achieved using spectra obtained with 
detectors placed far from the target (therefore subtend- 
ing small acceptance angles). Acceptance corrections are 
made using a calculation that combines Goldhaber's for- 
mulation of fragment angular distributions [17| with well- 
known Coulomb multiple scattering distributions, which 
are only important when high-Z targets are used (Sn, 
Pb) . The acceptance correction method has been shown 
to work well in our earlier published data [8l-[l0|. based 
on the fact that cross sections for heavier fragments mea- 



sured at large acceptance can be accurately reproduced 
by the corrected small-acceptance results. The main un- 
certainty in the calculation is associated with ao, the pa- 
rameter in the Goldhaber model that controls the widths 
of the momentum distributions of outgoing fragments af- 
ter the collision. The light-fragment cross sections allow 
model tests at large AZ (i.e., smaller impact parame- 
ters), where few previous comparisons have been made. 
This is a more stringent test of models than is possi- 
ble using data dominated by peripheral collisions (lead- 
ing to small charge changes). However, the limitations 
of the experimental method, which does not account for 
non-leading light fragments, and the meaning of the re- 
ported cross sections arc not entirely obvious and must be 
considered when making comparisons to models. These 
points will be further elucidated below. 

There are some instances of overlap between the mea- 
surements presented here and those made by Webber et 
al. [3 list ], and we present comparisons where data sets 
are sufficiently similar. Webber et al. obtained data on 
hydrogen and carbon targets, while our chosen targets 
span the periodic chart to allow for study of target mass 
dependences in the cross sections. 

The data presented here are drawn from a series of 
fragmentation experiments performed between 1995 and 
2006. Analysis of these data continues with the specific 
goal of extracting light fragment cross sections. A paral- 
lel effort is in progress to validate and verify the nuclear 
physics models used in space radiation shielding applica- 
tions [III ■ The progress of model development over time 
is being tracked by placing the models under configura- 
tion control, with automated validation benchmarks to 
enable comparisons as models are improved. Validation 
metrics are focused on the specific applications of inter- 
est, and have been developed to enable comparisons of 
fragmentation models to the relatively sparse experimen- 
tal database. The data shown here improve and extend 
that database. 



II. EXPERIMENTS 

The cross sections have been obtained from several sep- 
arate experiments which shared a common design. The 
16 data at 600 and 1000 MeV/nucleon were obtained 
at the NASA Space Radiation Laboratory (NSRL) at 
the Brookhaven National Laboratory. All other exper- 
iments were performed at the Heavy Ion Medical Ac- 
celerator in Chiba (HIMAC) at the Japanese National 
Institute of Radiological Sciences. In all cases, we iden- 
tify particles using deposited energy (AE) signals from 
small-area silicon detectors centered on and normal to the 
beam axis. Detectors are positioned just upstream of the 
target position so that event samples can be limited to 
those with one and only one well-identified primary beam 
ion present incident on the target, with a position close 
to the nominal beam axis. Other detectors are placed 
downstream of the target, at various distances with re- 
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spect to the target-center position so that they subtend 
different acceptance angles and measure different spec- 
tra. The large acceptance detectors, typically placed so 
as to subtend forward cones with half- angles between 5° 
and 10°, have 100% acceptance for surviving primaries 
and all fragments whose species can be identified, which 
generally extends as far as charges Zf lag > Z pT i maxy /2. 
For lower A£"s, there appears to be poor resolution, but 
this is in fact caused by the many possible combinations 
of light fragments, which results in overlapping AE dis- 
tributions. Detectors placed downstream so as to have 
small acceptance, on the order of 1° to 2°, are hit by 
many fewer fragments, and produce spectra in which all 
fragment species can be resolved. Given that the detec- 
tors are unsegmented, even at small acceptance there is 
some unavoidable ambiguity in interpretation of some of 
the fragment peaks, since certain combinations of light 
fragments can be indistinguishable from a single frag- 
ment of a heavier species (e.g., two He fragments with 
a particular velocity, detected in coincidence, produce a 
signal in the detector very close in amplitude to that of 
a single Li fragment at the same velocity). 

Figure 1 is a schematic drawing of the arrangement 
of detectors on the beamline for the 400 MeV/nucleon 
14 N experiment. It is representative of the configura- 
tions used for all experiments. A detailed discussion of 
the experimental setup can be found in Zeitlin et al. Q . 
The acceptance angles, defined as the half-angle of the 
forward cone extending from the target center to the ra- 
dius of the detector, are as indicated in the figure. The 
detectors were arranged in pairs to facilitate the data 
analysis, which depends on correlations between neigh- 
boring detectors. 



collection in one of the detectors. For analysis of small- 
acceptance detector data, we selected those events in 
which the charge of the most forward-produced fragment 
could not be determined by the large-acceptance detec- 
tors, plus events corresponding to the two or three light- 
est fragment species that can be resolved at large accep- 
tance. This provides overlap in the ranges of Z measured 
in the two acceptances, which in turn allows us in a sub- 
sequent step to tune the acceptance model so that the 
cross sections match as closely as possible for the frag- 
ment species measured at both acceptances. 



A. Analysis Using Large- Acceptance Detectors 

There is a degree of subjectivity involved in drawing 
the cut contours in the scatter plots. The contour that 
defines the surviving primary ions is particularly impor- 
tant, since it directly affects the charge-changing cross 
section <j cc . Since the fragment production cross sec- 
tions are effectively normalized to a cc , this cut strongly 
affects all measurements for a given data set. The 
contours define the event selection efficiencies, that is, 
N pass {Z) = N true (Z)e(Z), where N pass (Z) is the num- 
ber of events of a particular charge Z within the contour, 
and Ntrue(Z) is the "real" number of such particles. Our 
goal is to draw the contours so that all values of e(Z) are 
equal. However, there is no way to do this with perfect 
precision, and as a result the drawing of the contours is 
a source of systematic error that must be accounted for. 



B. Analysis Using Small- Acceptance Detectors 



III. EVENT SELECTION AND PARTICLE 
IDENTIFICATION 

With the beam energies and targets used in these ex- 
periments, primary ions that survive traversal of the tar- 
get and projectile fragments generally have velocities that 
fall in a narrow range. This results in AE spectra with 
peaks that are well separated from one another, with AE 
oc Z 2 . Both the number of visible fragment peaks in a 
particular spectrum, and the correspondence of the inte- 
grated counts in those peaks to cross sections, depend on 
the several factors described below. 

As in earlier work (see references @ - [13]), the CERN 
library program PAW [2(| was used to analyze the data. 
For every run, an initial event sample was selected by 
requiring that one and only one primary beam ion was 
seen in the detectors upstream of the target. For each 
detector pair downstream of the target, a scatter plot 
of AE's was made, and a cut contour (in some cases, 
multiple contours) was drawn to select those events hav- 
ing correlated pulse heights. These cuts remove events 
in which fragmentation occurred within either detector, 
along with events in which there was only partial charge 



For each run, we chose a subset of the events analyzed 
in the large-acceptance data for further analysis using 
the small acceptance detectors. The subsamples consist 
mostly of events in which the AE in the large-acceptance 
detectors is in the unresolved portion of the spectrum. 
The remainder of the events chosen are those in which 
the charge as determined by the large-acceptance detec- 
tors is at the low end of what can be resolved. Events 
with well-correlated signals in the downstream detector 
pair are selected by drawing another cut contour in the 
appropriate scatter plot. 

A typical large-acceptance fragment charge spectrum 
is shown in Figure 2a for the 400 MeV/nucleon 14 N data, 
with a small acceptance charge spectrum shown in Figure 
2b. The comparison illuminates some basic physics. In 
the large-acceptance data, only two fragment species, B 
and C, can be resolved at large acceptance. In contrast, 
at small acceptance, considerable structure is visible in 
the charge histogram. The scale in this plot was de- 
termined simply by scaling the square root of deposited 
energy so that the peak for C fragments (not shown) 
is centered at a charge of 6.0. As in previously-reported 
data sets, non-integer peaks are seen in addition to peaks 
for charges 1 through 5. There is a peak near 4.5 due to 
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FIG. 1. Schematic diagram of the beamline configuration for the 400 MeV/nucleon 14 N experiment. Spacing between detectors 
was 2 cm unless otherwise noted. Configurations for other experiments were similar. 



detection of Be and He fragments in coincidence. There 
is also a broad peak centered near 3.5, due to the de- 
tection of three helium nuclei in coincidence. Since the 
beam ion was 14 N, charge conservation allows for a pro- 
jectile fragment of charge 1 ( 1 H or 2 H) to also be present 
in these events. (The presence or absence of a charge 1 
fragment contributes to the width of the Z«3.5 peak.) 
The peak in the vicinity of charge 3 appears to be split, 
with a relatively sharp peak centered very close to 3.0 
and a less-defined peak near 2.8. The former is likely 
due to Li fragments, and the latter to pairs of He frag- 
ments detected in coincidence, some of which come from 
decay of any 8 Be fragments. In most of our other data 
sets, the two peaks near charge 3 are not resolved, and 
though the separation here is not large, there do seem to 
be two distinct peak regions. Fitting two Gaussian dis- 
tributions to this region yields a x 2 of 12.1 for 10 degrees 
of freedom, with one peak found at Z = 2.76 ± 0.013 and 
the other at 2.96 ± 0.013 [Hj]. Between charge 1.8 and 
about 2.6, there are three statistically-significant peaks 
at 1.94, 2.15, and 2.46, likely corresponding to He alone, 
and He in coincidence with one and two H fragments, 
respectively. In the vicinity of charge 1, peaks appear 
at 0.88 and 1.36, corresponding to one and two singly- 
charged relativistic particles. Another notable feature of 
this spectrum is the suppression of charge 4 fragments. 
As mentioned above, any 8 Be fragments created in the 
collisions instantaneously decay to two 4 He ions, which 
may then be detected in the Z » 3.5 or Z w 2.8 peaks, 
both of which are far more populated than the charge 4 
region, which extends from about 3.8 to 4.7. 

Additional complications arising in the interpretation 
of small-acceptance spectra are discussed elsewhere [22j ■ 
It was noted in that article that, at least for one sim- 
ulated data set (650 MeV/nucleon 40 Ar), some 80% of 
He fragments were produced in association with heavier 
fragments. More generally, we can say that significant 
percentages of the lightest fragments are produced as 
non-leading fragments, and are difficult to account for. 
We note here that the light-fragment results presented 
below are quite sensitive to the acceptance angles in the 
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FIG. 2. Charge spectra at large (upper, Fig. 2a) and small 
(lower, 2b) acceptances, for 400 MeV/nucleon 14 N beam on a 
carbon target. 



different experiments, and that large correction factors 
are applied. The net result is that relative errors on the 
cross sections for the lightest fragments are large com- 
pared to those obtained for heavier fragments. 

With the 14 N beam, interpretation of the multiple- 
fragment peaks is relatively straightforward. Charge con- 
servation dictates that (unlike in data sets with heavier 
beams) the charge 4 peak cannot be due to the detec- 
tion of four He fragments in coincidence. When counting 
events by leading fragment, as is done in the cross sec- 
tion analysis, the events populating the peak near 3.5 are 
counted as He, and the only ambiguity is in assigning a 
charge to the events in the 2.6-3.5 region, as they may be 
either Li fragments or pairs of He fragments. For simplic- 
ity, we group them together in the cross section tables, 
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but separately we will show estimates of the He:Li ratios. 



C. Experiment Acceptances 

The detector configuration changed between experi- 
ments, which were conducted over a period of several 
years at different accelerator facilities. Table I shows the 
detector acceptance angles for the data sets presented 
here. The acceptance angles are the half-angles of the 
forward cones as measured from the exact center of the 
target. 

The variation in large acceptance angles is not signif- 
icant. In all cases, acceptance of fragments with Z > 
■Zbcam/2 is calculated to be at least 99%. The variations 
in small acceptance angles are more significant, as the 
measured spectra are found to be quite sensitive to both 
angle and beam energy. 



D. Correction Factors 

The raw counts of events by species must be corrected 
for various effects before cross sections can be computed. 
The correction factors have been described in detail pre- 
viously [zlGjJ. The corrections are summarized in Table 
II. All have been applied to the cross sections presented 
below. 

The magnitude of a given correction depends on the 
depth of the target (due to secondary and higher-order 
interactions in the target) and the configuration of the 
silicon stack in a particular run. Corrections are small- 
est for the case of large-acceptance detectors close to a 
thin target. The correction factors have associated un- 
certainties that are taken into account when estimating 
systematic errors. Of the three effects listed in Table II, 
the corrections for multiple interactions in the target are 
typically the largest, particularly for the lower-mass tar- 
gets (CH2 and C) where even modest depths on the order 
of 3 g cm -2 cause « 20% of the beam ions to fragment, 
leading to « 10% corrections for the heaviest fragments. 

Projectile fragments generally receive fairly small 
transverse momenta in the collisions, so that their an- 
gular distributions are strongly forward-peaked. Multi- 
ple scattering angles are small at the energies considered 
here, and can be ignored in determining large-acceptance 
cross sections. Because, as mentioned above, the large- 
acceptance data also require the smallest corrections for 
losses in the detectors and intervening materials, we use 
them to obtain the charge-changing cross sections (cr cc ) 
and fragment cross sections for all resolvable species. As 
can be seen in Fig. 2 above, where only peaks for B and C 
fragments are clear, this can be as little as two fragment 
species. 



IV. SYSTEMATIC ERRORS 

The statistical errors are generally small in these exper- 
iments, but the systematic uncertainties contain several 
contributions and typically dominate the total. When 
sufficient beam time is available, we take at least two runs 
with the same beam ion/energy /target combinations and 
vary the depth of the target. Cross sections obtained at 
different depths of the same material must, after cor- 
rections, be equal. The variations in the cross sections 
obtained this way are a good measure of the overall sys- 
tematic uncertainties. When we do not have multiple 
data sets to combine, we associate a conservatively-large 
error with our data selection cuts and propagate the un- 
certainties into all cross sections, as described in the next 
section. 

As a practical matter, one can only obtain reasonable 
fragment statistics by using targets whose depths rep- 
resent at least a few percent of an interaction length. 
With high-Z targets like Sn and Pb, ionization energy 
losses per unit interaction length are large compared to 
low-Z materials. In order to keep the beam energy ap- 
proximately constant throughout the depth, high-Z tar- 
gets must therefore be kept thin, yielding poorer frag- 
ment statistics and larger relative systematic errors on 
the cross sections compared to lower-mass targets such 
as C and Al, for reasons that will become apparent in the 
following. 



A. Uncertainties in the Charge-Changing Cross 
Sections 

The definition of the cut contour that defines the sur- 
viving primaries and heaviest fragments is the single 
largest source of systematic error in these measurements. 
The most difficult contours to draw are those for runs 
with either no target or a very thin target, because in 
these cases the tail of the primary distribution on the 
low side can be difficult or impossible to distinguish from 
AZ = 1 events. Even with thicker targets, there is al- 
ways at least a small number of ambiguous events that 
fall between the clusters of primary ions and the highest- 
Z fragments. To account for the cut contour uncertainty, 
a systematic error Sf is assigned to the fraction / of sur- 
viving primaries. The magnitude of Sf is determined by 
repeatedly drawing the contours and examining the re- 
sults. Typical variations in / are 0.005 or smaller. In 
order to get a better sense of the magnitude of this un- 
certainty as it propagates into the charge-changing cross 
section, er cc , consider that the cross section scales with 
ln(f). For thin targets, it is approximately true that 
a cc oc (1 - /), so that 6a cc /a cc oc 5f/(l - f). Since Sf is 
found to be more or less constant for a given experiment, 
5a cc is largest when / approaches 1 , corresponding to the 
thin-target case. 

The preceding argument would appear to favor the use 
of thicker targets. However, that is only true to the ex- 
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TABLE I. Acceptance angles for all data sets. 



Beam Ion 


Energy at Extraction 
(MeV/nucleon) 


Large Acceptance Angle 
(degrees) 


Small Acceptance Angle(s) 
(degrees) 


i4 N 


290 


5.7° 


1.7° 


14 N 


400 


9.8° 


3.8° 


16Q 


290 


5.7° 


1.7° 





400 


6.7° 


2.5° 


16Q 


600 


7.6° 


1.4° 


16Q 


1000 


7.6° 


1.4° 


20 Ne 


290 


4.8° 


1.6° 


20 Ne 


400 


7.0° 


1.8° 


20 Ne 


600 


7.0° 


2.5°, 1.7° 


24 Mg 


400 


9.5° 


2.0° 



TABLE II. List of correction factors applied to charge spectra. 



Physical Effect 



Effect on Spectra 



Correction Estimation Method 



Charge-changing interactions in 
air gaps, dead layers, etc. 
Multiple interactions in target 
Charge-changing interactions in 
the detector stack 



Loss of primary ions and 
increased numbers of fragments 
Shifts of fragment distributions 
Z-dependent detection efficiency 



Target-out data 

Monte Carlo simulation [23] 
Geometric cross section model 



tent that cross sections can be approximated as energy 
independent. As target depth increases, so too does the 
range of energies at which the primaries can interact. 
This can obscure subtleties in the energy dependences of 
either a cc or the fragment production cross sections. 



The model used to correct for losses due to nuclear in- 
teractions in the detectors is estimated to contribute on 
the order of ± 1% relative error to the systematic uncer- 
tainty in cr cc , and the accuracy of the target areal den- 
sity measurements is also estimated at ± 1%. For thin 
targets, the uncertainty associated with the cut contour 
dominates the quadrature sum, but these smaller contri- 
butions can be important for thicker targets. Determi- 
nation of the systematic error on a given a cc is made by 
combining the results from multiple targets. In combin- 
ing data sets, the weighted average and a \ 2 are com- 
puted. The error on each measurement is initially set 
equal to that arising from the definition of the primary 
selection cut contour. If this initial \ 2 is found to be 
greater than 1.0 per degree of freedom (the number of 
data sets minus one) , then an additional systematic error 
is added in quadrature to the starting errors of each in- 
dividual measurement, and incremented upward in steps 
of 0.1% relative error until we achieve the desired % 2 re- 
sult (< 1 per degree of freedom). In practice, these extra 
contributions are often not required since the estimates 
of the cut contour uncertainties are conservatively large 
and the initial calculation of \ 2 typically yields a value 
< 1 per degree of freedom. 



B. Uncertainties on Fragment Cross Sections at 
Large Acceptance 

For each data set, the uncertainty on the charge- 
changing cross section is propagated into the fragment 
cross sections, and added in quadrature to the statistical 
errors. Statistical errors are much more significant for 
fragments than for primaries. When data sets are com- 
bined, we again allow for the addition in quadrature of 
additional systematic uncertainties sufficient to bring the 
total \ 2 f° r combining data (summed over all fragment 
species) to less than 1.0 per degree of freedom. We find 
it is common that non-zero addition uncertainties are re- 
quired at this step, in contrast to what is seen when com- 
bining data sets to obtain charge- changing cross sections. 
We believe the major contributions that these ad hoc ad- 
ditions are covering are the uncertainties associated with 
ambiguities in (1) the cut contour that defines the sample 
of the heaviest fragments, and (2) the counting of events 
in the "valleys" between fragment peaks. 



C. Additional Uncertainties on Fragment Cross 
Sections at Small Acceptance 

The uncertainties defined above contribute to the light- 
fragment cross sections. In addition, the acceptance cor- 
rections that are made to account for the fragment an- 
gular spreads contribute to the overall uncertainty. Pre- 
viously, the relative errors associated with these correc- 
tions were found to be about ± 5-6% 10]. This accounts 
for reasonable variations in the size of the beam and its 
divergence, the model parameter ao, which controls the 
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widths of the momentum distributions, and the exact ac- 
tive areas of the detectors. Here, a ± 6% uncertainty is 
added in quadrature with the other uncertainties. 

V. CHARGE-CHANGING CROSS SECTIONS 

Table III shows the charge-changing cross section re- 
sults. The same results are shown in Figure 3, but with- 
out the hydrogen-target cross sections in order to keep 
the y-axis readable, and with scaling factors applied as in- 
dicated in the figure legend, again for reasons of readabil- 
ity. NUCFRG2 and PHITS predictions are also shown 
in Table III. The model comparisons are discussed in the 
following section. 

In general, the charge-changing cross sections shown 
here for carbon and heavier targets follow the same trends 
reported in our previous work with other beam ions in 
the same energy range: there is little or no energy de- 
pendence, and a simple geometrical model of overlap- 
ping spheres with a transparency term reproduces the 
data well. The only apparent oddity in the results in 
Table III is the fact that the two hydrogen-target cross 
sections for 14 N beams at 290 and 400 MeV/nucleon are 
larger than those obtained with 16 O beams at the same 
energies. However, that refers only to the central values; 
taking into account the uncertainties, the data are al- 
most equally consistent with the hypothesis that the 16 O 
cross sections are larger, as one would expect from sim- 
ple geometry. The only other unexpected trend in these 
data is the slight decrease of the 20 Ne cross sections for 
aluminum and heavier targets at 600 MeV/nucleon com- 
pared to the 400 MeV/nucleon data. The effect is only 
slightly larger than the combined one-sigma uncertainties 
and we do not believe it has any physical significance. 



A. Comparison to Previous Measurements 

The most relevant published data to which we com- 
pare are drawn from a large number of charge-changing 
cross sections [18f and fragment production cross sections 
14] published by Webber et al. We have made sev- 
eral previous comparisons to these data; in many cases, 
the differences are larger than the uncertainties. For in- 
stance, in Ref. Q, we found discrepancies on the order of 
5% between our data and Webber's when both carbon- 
target and hydrogen-target data were considered. Larger 
discrepancies were seen in several instances where frag- 
ment production cross sections were compared. The sit- 
uation is similar with the present data sets. Table IV 
shows results for seven of the ten data sets analyzed here; 
the seven were selected for similarities in beam energies. 
Agreement between the hydrogen-target cross sections 
is generally good, but it is not good for several of the 
carbon-target results. Specifically, the cross sections for 
oxygen beams on carbon targets do not agree well, and 
the results for 20 Ne on carbon at 400 MeV/nucleon dif- 



fer by about 8%. The discrepancies are all well beyond 
the stated uncertainties from either group. An additional 
independent measurement might be useful in these cases. 

Since both groups used polyethylene and carbon tar- 
gets to obtain the hydrogen-target results, and those are 
in better agreement, it must also be the case that the 
polyethylene-target results have a level of disagreement 
similar to that seen for the carbon targets. The agree- 
ment for hydrogen is to some extent a product of the 
cancellation of errors in the subtraction. 

In the results shown here, the level of agreement be- 
tween our experimental results and those of Webber et 
al. is not especially good. If errors were correctly es- 
timated in all experiments, we would expect to see very 
few values of x 2 as large as those seen in Table IV (i.e., > 
5 for one degree of freedom). This comment pertains to 
several previously-published data points as well. In some 
instances, there are additional data from other groups 
that can be used for additional comparisons. Our charge- 
changing cross sections are generally in good agreement 
with those measured using plastic nuclear track detectors 
00. 



B. Charge-changing Cross Sections Compared to 
Model Predictions 

A discussion of the physics content of the various mod- 
els used here is given below in Section IX. A. 

As in earlier work, we compare the charge-changing 
cross section data to several models, including a simple 
formula given by a cc — lOirrQ^A]^ 3 +Ap^ 3 — 6) 2 , allowing 
b (dimensionless) and ro (in units of fm) to vary. Where 
A t and A p are the atomic mass numbers of the target 
and projectile, respectively. For each charge-changing 
cross section and each model, we calculated the ratio of 
the predicted cross section to the measured. We then 
averaged the results over all beam ion/energy combina- 
tions for a particular target species. We then examined 
the aggregate results, that is, the results for all targets 
combined. Tuning the parameters to values of b = 0.938 
and ro = 1.375 fm yields results that, using this aver- 
aging method, agree with the data to better than 0.3%. 
(For the 600 MeV/nucleon 20 Ne data reported earlier, we 
found good agreement with the data for b = 0.86 and ro 
= 1.34 fm.) The near-perfect average agreement is not 
particularly significant - the standard deviation relative 
to the mean is a better measure of typical accuracy over 
multiple beam ions, energies, and targets. 

The averaged results are shown in Table V, for this 
simple form and for NUCFRG2 and PHITS as well. 
Values below 1.0 indicate that the predicted cross sec- 
tions are smaller than the measured values. The simple, 
energy-independent model and PHITS yield the same 
value (0.041) for the figure of merit, whereas NUCFRG2 
gives a somewhat higher value (0.050), indicating more 
variance between the data and NUCFRG2 than for the 
other models. If we exclude the cross sections obtained 
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FIG. 3. Scaled charge-changing cross sections for all beams reported here, for targets from carbon to lead. 



TABLE III. Measured charge-changing cross sections and uncertainties, along with NUCFRG2 and PHITS predictions. Listed 
energies are for the extracted beams. 



Beam Ion, Energy 


H target 


C target 


Al target 


Cu target 


Sn target 


Pb target 


14 N, 290 


229 ± 11 


750 ± 14 


1190 ± 31 


1794 ± 64 


2303 ± 140 


3092 ± 197 


NUCFRG2 


223 


853 


1198 


1741 


2353 


3046 


PHITS 


231 


782 


1149 


1767 


2314 


3016 


14 N, 400 


236 ± 11 


795 ± 19 


1170 ± 36 


1784 ± 45 


2394 ± 78 


3193 ± 88 


NUCFRG2 


232 


866 


1215 


1764 


2384 


3090 


PHITS 


246 


776 


1136 


1718 


2227 


2985 


ib O, 290 


219 ± 13 


863 ± 20 


1219 ± 41 


1798 ± 60 


2369 ± 74 


3125 ± 118 


NUCFRG2 


226 


923 


1287 


1857 


2495 


3217 


PHITS 


223 


832 


1220 


1859 


2436 


3231 


ib O, 400 


220 ± 17 


842 ± 22 


1186 ± 27 


1772 ± 51 


2512 ± 72 


3121 ± 74 


NUCFRG2 


237 


937 


1304 


1881 


2528 


3264 


PHITS 


245 


827 


1208 


1840 


2434 


3217 


ib O, 600 


264 ± 17 


902 ± 28 


1206 ± 29 


1892 ± 40 


2524 ± 69 


3366 ± 112 


NUCFRG2 


259 


976 


1351 


1940 


2603 


3359 


PHITS 


269 


856 


1211 


1818 


2404 


3224 


iB o, 1000 


276 ± 16 


856 ± 26 


1260 ± 19 


1864 ± 40 


2518 ± 69 


3307 ± 112 


NUCFRG2 


285 


1024 


1408 


2011 


2692 


3479 


PHITS 


280 


937 


1211 


1771 


2312 


3116 


^ u Ne, 290 


272 ± 16 


1050 ± 21 


1445 ± 45 


2043 ± 99 


2807 ± 142 


3556 ± 260 


NUCFRG2 


287 


1062 


1462 


2083 


2772 


3548 


PHITS 


281 


998 


1420 


2107 


2776 


3543 


^ u Ne, 400 


311 ± 15 


1034 ± 21 


1438 ± 26 


2140 ± 47 


2764 ± 99 


3555 ± 129 


NUCFRG2 


298 


1078 


1481 


2109 


2810 


3603 


PHITS 


305 


983 


1416 


2092 


2723 


3538 


2U Ne, 600 


319 ± 13 


986 ± 14 


1349 ± 30 


1993 ± 35 


2572 ± 75 


3407± 73 


NUCFRG2 


321 


1119 


1532 


2173 


2892 


3711 


PHITS 


326 


1018 


1422 


2098 


2758 


3580 


a Mg, 400 


328 ± 18 


1028 ± 18 


1480 ± 38 


2244 ± 79 


2794 ± 142 


3727 ± 239 


NUCFRG2 


315 


1124 


1547 


2194 


2935 


3761 


PHITS 


313 


998 


1428 


2117 


2808 


3670 
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with 290 MeV/nucleon beams, the simple model yields 
an even better figure of merit, 0.035. No similar improve- 
ment is seen if the 290 MeV/nucleon data are excluded 
from the PHITS and NUCFRG2 comparisons. Thus, the 
simple energy independent model can, at least for this 
very limited range of beam ions and energies, be tuned 
to a high degree of accuracy. It seems likely that other 
parameter values can be found that would give better 
agreement with other data sets (higher-energy and/or 
higher-mass beams); thus it may be reasonable to treat 
one or both of these parameters as energy dependent to 
obtain an empirical fit across a greater range of data. 

Some trends are apparent in the averages by target. 
For NUCFRG2, the ratios of predicted to measured cross 
sections are significantly greater than 1 for carbon and 
aluminum targets, with values closer to 1.0 found for Cu, 
Sn, and Pb targets. NUCFRG2 also does well, on aver- 
age, for hydrogen targets. The lack of agreement with 
the data for C and Al targets is potentially important, 
since the beam ions studied here are relatively abundant 
among GCR heavy ions, and both C and Al are impor- 
tant shielding materials in space. For PHITS, the aver- 
ages are highly consistent for C, Al, and Cu, then de- 
crease slightly for Sn and Cu. 

The comparatively small spread in the ratios found 
for the simple model does not hold when this formula 
with these same parameter values is applied to other 
beam ion/target combinations. When we average over 
our published data for beams from 290 MeV/nucleon 12 C 
to 1000 MeV/nucleon 56 Fe, the average remains close to 
1.0 (0.994), but the standard deviation increases to about 
6.4%, and some discrepancies appear that approach a 
10% level of disagreement. Many of these data points 
were obtained at higher energies than the data presented 
here. 

In the following we will compare fragment production 
cross sections not only with PHITS and NUCFRG2, as 
above, but also with the predictions of the Los Alamos 
version of the Quark Gluon String Model (LAQGSM) 
[H In the output generated by the LAQGSM code, 
charge-changing cross sections are not provided. Rather, 
cross sections are computed as elastic or inelastic. This 
does not allow for comparisons to the data above, since 
the inelastic and charge-changing cross sections do not 
correspond to one another. The inelastic cross section 
contains the charge-changing cross sections plus contri- 
butions from neutron-stripping reactions, and is therefore 
larger than the charge-changing cross section alone. The 
clastic cross section cannot be measured in our experi- 
ments. Therefore, meaningful comparisons of the charge- 
changing cross section data to LAQGSM predictions are 
not possible given the current state of the output from the 
code. But as we will show below, valuable comparisons 
can be made for fragment production cross sections. 



VI. FRAGMENT PRODUCTION CROSS 
SECTIONS AT LARGE ACCEPTANCE 

Fragment production cross sections are presented in 
two parts. We first discuss those obtained at large ac- 
ceptance, where it is typically possible only to measure 
fragment species with (approximately) Z[ rag > Z\> ca _ nl /2. 
These data are readily compared to model predictions, 
and in some instances, we are also able to c omp are to pre- 
viously published data from Webber et al. [Tj] for hydro- 
gen and carbon targets. In the second part of the discus- 
sion, we present cross sections for lighter fragments mea- 
sured with the small-acceptance detectors. These cross 
sections are (with a few exceptions) corrected for accep- 
tance losses using a previously-described method 0, II] 
based on Goldhaber's model of nucleon momentum dis- 
tribution within the nucleus [171 . subsequently modified 
by Tripathi and Townsend [27] . Complications in inter- 
preting these light-fragment cross sections are discussed 
below. 

Direct comparisons of the results obtained at large ac- 
ceptance can be made to a few previous measurements, 
particularly those of Webber et al. [HI • We return to this 
point below in Section IX. E. No measurements compara- 
ble to our small-acceptance results are in the literature 
to our knowledge. 



A. Nitrogen Beams 

For the 290 MeV/nucleon 14 N data, we are able to re- 
solve charge 4 peaks, whereas in the 400 MeV/nucleon 
14 N spectra we can only see fragment peaks for charges 
5 and 6. This is probably due to the greater forward- 
focusing of non-leading fragments at the higher energy, 
resulting in a higher average detected multiplicity. Table 
VI shows the results for nitrogen beams at two energies. 
The charge 4 cross sections for the 290 MeV/nucleon 
beam are presented below, with the small-acceptance re- 
sults. The results are quite similar at these two ener- 
gies. The AZ = 1 cross sections appear systematically 
smaller by 8-10% at 400 MeV/nucleon for all but the H 
and Pb targets. (For the Pb target, the large error on 
the 290 MeV/nucleon data point precludes drawing any 
such conclusion.) This is qualitatively consistent with 
trends we have observed in other data sets, in which the 
cross sections for the smallest charge changes decrease 
slightly with energy while those for larger charge changes 
are more constant or even increase. For boron fragments 
(charge 5) no significant differences are seen between the 
two energies. The boron production cross sections are 
roughly a factor of two lower than those for carbon pro- 
duction. For charge 4 fragments, the pattern seen for 
charge 6 fragments repeats; that is, for C, Al, Cu, and 
Sn targets, the 400 MeV/nucleon cross sections are sys- 
tematically smaller by 10-25% than those measured with 
the 290 MeV/nucleon beam. For both charge 4 and 6 
fragments, the Pb-target uncertainties are too large to 
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TABLE IV. Charge-changing cross sections for hydrogen and carbon targets as reported here and as measured by Webber ct 
al. Listed beam energies are at the target centers. The x 2 values are all for one degree of freedom. 



Beam Ion 


Energy 


H target 


H target 


X 2 for 


C target 


C target 


X 2 for 




(MeV/nucleon) 


This work 


Ref. [18] 


agreement 


This work 


Ref. [18] 


agreement 


N 


375, 516 


236 ± 12 


227 ± 6 


0.5 


795 ± 19 


796 ± 8 


0.00 





375, 441 


224 ± 17 


232 ± 6 


0.2 


872 ± 22 


794 ± 8 


11.1 





578, 491 


264 ± 17 


247 ± 6 


0.9 


902 ± 28 


823 ± 8 


7.4 





980, 903 


276 ± 16 


248 ± 6 


2.7 


856 ± 26 


813 ± 8 


2.5 


Ne 


375, 468 


311 ± 15 


298 ± 7 


0.6 


1034 ± 21 


951 ± 10 


12.7 


Ne 


572, 599 


311 ± 14 


319 ± 8 


0.2 


984 ± 13 


980 ± 10 


0.1 


Mg 


358, 309 


328 ± 18 


310 ± 8 


0.8 


1028 ± 19 


1097 ± 11 


9.9 



TABLE V. Averaged ratios of predicted charge-changing cross sections to measured cross sections. There are two entries in each 
cell, the first being the averaged ratio and the second being the standard deviation. The simple model is energy-independent, 
while NUCFRG2 and PHITS contain slight energy dependence, which is weak except for hydrogen-target cross sections. In the 
bottom row, ratios have been re-computed for the simple model excluding data from the 290 MeV/nucleon beams. 



Model 


H target 


C target 


Al target 


Cu target 


Sn target 


Pb target 


Combined 


NUCFRG2 
PHITS 
Simple 
Simple* 


1.006, 0.041 
1.021, 0.041 

n/a 

n/a 


1.097, 0.052 
0.991, 0.049 
1.019, 0.055 
1.015, 0.041 


1.066, 0.048 
0.991, 0.029 
0.995, 0.037 
1.000, 0.036 


1.023, 0.043 
0.994, 0.041 
0.980, 0.036 
0.975, 0.040 


1.036, 0.041 
0.985, 0.046 
1.003, 0.036 
0.997, 0.022 


1.019, 0.036 
0.982, 0.034 
1.010, 0.033 
1.010, 0.030 


1.041, 0.050 
0.994, 0.041 
1.000, 0.041 
1.000, 0.035 



say whether or not they are consistent with the trend 
observed for C through Sn targets. 

In the last two rows of Table VI, we show the cross 
section for events in which the charge of the leading par- 
ticle is less than 4. Using the large acceptance detectors, 
there is no more information to be gleaned from these 
events. However, it is clear that these cross sections in- 
crease monotonically with target mass as a share of the 
charge changing cross section. These cross sections gen- 
erally increase with increasing beam energy. The trend is 
not statistically significant in the 14 N data owing to the 
small difference in the two beam energies in these data, 
but it is seen clearly in other data sets. 



B. Oxygen Beams 

For the four 16 O beams studied, fragment charges from 
5 to 7 can be resolved at large acceptance. The results are 
shown in Table VII. As in the 14 N data above, the AZ = 
1 cross sections (charge 7) tend to decrease with increas- 
ing beam energy for all targets except H and Pb. This 
is also predominantly the case for AZ = 2 and AZ = 3, 
though there are a few exceptions to the general trend. 
Since the charge-changing cross sections are, except for 
hydrogen targets, roughly constant, the cross sections for 
the category Z < 5 increase with increasing beam energy. 
This is a modest effect, on the order of 10 to 15% for car- 
bon and heavier targets. The hydrogen target data show 
somewhat different behavior: the charge-changing cross 
section increases with energy, and the cross section in the 
Z < 5 category also increases as energy increases, while 
the cross sections for AZ = 1, 2, and 3 are approximately 



independent of energy in this range. 

The odd-Z even-Z effect, discussed in detail elsewhere 
[loj . can be seen here: the production cross sections for 
AZ = 2 are higher than the corresponding cross sections 
for AZ = 1 in all cases for H, C, Al, and Cu targets, 
but not for Sn and Pb targets. A similar trend was seen 
in silicon-beam data Q. However, with heavier beams 
( 40 Ar, 48 Ti, and 56 Fe), the opposite is true. As was the 
case for the 14 N beams, the boron production cross sec- 
tions are 50 to 60% as large as those for carbon produc- 
tion. 



C. Neon Beams 

Three 20 Ne beam energies were studied. Fragment 
charges from 5 to 9 can be resolved at large acceptance. 
The results are shown in Table VIII. The same four 
trends discussed above for N and O beams can be seen 
here: even-Z fragment production cross sections are en- 
hanced compared to those for odd Z's; for a given frag- 
ment species and target, cross sections for the species 
that can be resolved decrease with increasing energy for 
carbon and heavier targets; boron production cross sec- 
tions are again about a factor of 2 smaller than those for 
carbon production; and the cross section for the unre- 
solved category increases with increasing beam energy. 

We previously published results for 600 MeV/nucleon 
20 Ne on the same target materials 0. Those cross 
sections have been combined with an additional 600 
MeV/nucleon data set obtained subsequently. Our meth- 
ods for combining data sets and estimating systematic er- 
rors have changed in the interim; the results shown here 



11 



TABLE VI. Fragment production cross section for N beams on elemental targets. Beam energies at extraction are shown, in 
units of MeV/nucleon. All cross sections are in millibarns. 



■Zfrag 


■E'bcara 


H 


C 


Al 


Cu 


Sn 


Pb 


6 


290 


82 ± 8 


162 ± 3 


221 ± 7 


283 ± 11 


319 ± 22 


389 ± 28 


6 


400 


83 ± 4 


150 ± 4 


195 ± 8 


260 ± 8 


312 ± 13 


415 ± 15 


5 


290 


29 ± 3 


72 ± 2 


98 ± 3 


132 ± 6 


152 ± 11 


176 ± 13 


5 


400 


27 ± 2 


75 ± 2 


96 ± 4 


127 ± 4 


154 ± 7 


176 ± 7 


4 


290 


14 ± 2 


48 ± 2 


67 ± 3 


97 ± 5 


115 ± 9 


133 ± 11 


4 


400 


15 ± 2 


43 ± 2 


59 ± 4 


81 ± 5 


92 ± 7 


130 ± 9 


<4 


290 


104 ± 11 


468 ± 14 


804 ± 31 


1282 ± 64 


1717 ± 140 


2394 ± 197 


<4 


400 


111 ± 12 


526 ± 19 


820 ± 36 


1316 ± 45 


1836 ± 78 


2472 ± 88 



TABLE VII. Fragment production cross section for O beams on elemental targets. Beam energies at extraction are shown, 
in units of MeV/nucleon. All cross sections are in millibarns. 



■Zfrag 


■E'bcam 


H 


C 


Al 


Cu 


Sn 


Pb 


7 


290 


59 ± 3 


137 ± 4 


169 ± 7 


204 ± 8 


247 ± 10 


297 ± 15 


7 


400 


58 ± 4 


124 ± 4 


145 ± 5 


187 ± 8 


244 ± 11 


269 ± 10 


7 


600 


65 ± 4 


120 ± 5 


138 ± 4 


188 ± 6 


240 ± 10 


285 ± 14 


7 


1000 


59 ± 3 


105 ± 4 


131 ± 3 


167 ± 5 


210 ± 8 


270 ± 13 


6 


290 


61 ± 4 


148 ± 4 


181 ± 7 


225 ± 8 


258 ± 11 


305 ± 15 


6 


400 


60 ± 4 


131 ± 5 


159 ± 5 


209 ± 8 


241 ± 12 


257 ± 10 


6 


600 


71 ± 4 


131 ± 5 


154 ± 5 


206 ± 7 


248 ± 10 


279 ± 14 


6 


1000 


63 ± 3 


121 ± 4 


153 ± 4 


199 ± 6 


225 ± 8 


281 ± 13 


5 


290 


23 ± 2 


83 ± 2 


101 ± 4 


143 ± 5 


157 ± 13 


201 ± 11 


5 


400 


25 ± 2 


70 ± 3 


86 ± 3 


99 ± 5 


128 ± 7 


149 ± 6 


5 


600 


25 ± 2 


72 ± 3 


76 ± 3 


112 ± 4 


123 ± 6 


168 ± 9 


5 


1000 


29 ± 2 


56 ± 2 


72 ± 2 


91 ± 3 


114 ± 5 


130 ± 7 


<5 


290 


76 ± 13 


495 ± 20 


768 ± 41 


1226 ± 60 


1707 ± 74 


2322 ± 118 


<5 


400 


77 ± 17 


517 ± 22 


796 ± 27 


1277 ± 51 


1899 ± 72 


2446 ± 74 


<5 


600 


103 ± 13 


579 ± 28 


838 ± 29 


1386 ± 40 


1913 ± 69 


2634 ± 112 


<5 


1000 


125 ± 16 


574 ± 26 


904 ± 19 


1407 ± 40 


1969 ± 69 


2626 ± 112 



fully incorporate the newer methodology. 



D. 24 Mg at 400 MeV/nucleon 

Data were obtained with a 400 MeV/nucleon 24 Mg 
beam. Fragment charges from 6 to 11 can be resolved at 
large acceptance. These results are shown in Table IX. 
Also shown in Table IX are the cross sections for boron 
fragments (charge 5), which were obtained at small ac- 
ceptance; corrections have been applied. 



E. The Odd-Z Even-Z Effect 

For the 14 N and 16 beams discussed above, cross sec- 
tions for just three fragment species can be measured at 
large acceptance. This is not a sufficient number to al- 
low us to make any statements about the odd- Z /even- Z 
effect. However, with Ne and Mg beams, we can explore 
the magnitude of the effect in a quantitative way. Both 
beam ions have isospin T z = 0. As in previous work [Toj . 



we use the quantity V(Z) defined by Iancu et al. (l6| : 

V (Z f ) = 2a (Z f )/[a (Z f + l) + a(Z f - 1)] 

where Zf refers to fragments of charge Z. In Ref. (Toj . 
we combined results for C and Al targets. Here, exami- 
nation of the 20 Ne and 24 Mg data shows no statistically 
significant differences between C, Al, Cu, and Sn tar- 
gets, so all are combined. Lead (Pb) targets are excluded 
since results may be distorted by the contributions from 
electromagnetic dissociation, which produces increases in 
the AZ — 1 (and perhaps AZ = 2) cross sections. We 
combine the values of ViZf) obtained for all odd-Z frag- 
ments into a single weighted-average value, and similarly 
combine the results for all even-Z fragments to get that 
weighted average, and take the ratio of the two to obtain 
a single value for a given beam ion and energy. 

For 20 Ne beams, we can calculate V(Zf) only for 
charges 6, 7, and 8. We find overall ratios of 1.71 ± 
0.05, 1.73 ± 0.04, and 1.86 ± 0.04 for 290, 400, and 600 
MeV/nucleon beam energies, respectively. In Ref. [lCj . 
we used this same method for 28 Si beams at extracted 
energies of 290, 400, 600, 800, and 1200 MeV/nucleon, 
and found that the even:odd V[Zf) ratio for C and Al 
targets increased slightly as beam energy increased in 
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TABLE VIII. Fragment production cross section for Ne beams on elemental targets. Beam energies at extraction are shown, 
in units of MeV/nucleon. All cross sections are in millibarns. 



■^frag 


■^beam 


H 


C 


Al 


Cu 


Sn 


Pb 


9 


290 


49 ± 3 


109 ± 3 


130 ± 5 


175 ± 11 


229 ± 14 


308 ± 27 


9 


400 


47 ± 3 


102 ± 3 


120 ± 3 


160 ± 5 


195 ± 12 


236 ± 13 


9 


600 


53 ± 2 


84 ± 2 


106 ± 3 


144 ± 4 


177 ± 7 


235 ± 8 


8 


290 


78 ± 4 


163 ± 4 


188 ± 7 


230 ± 13 


281 ± 17 


304 ± 27 


8 


400 


81 ± 4 


151 ± 4 


178 ± 5 


228 ± 7 


262 ± 15 


318 ± 15 


8 


600 


75 ± 3 


134 ± 3 


158 ± 5 


205 ± 6 


237 ± 9 


296 ± 10 


7 


290 


54 ± 3 


128 ± 4 


155 ± 6 


191 ± 11 


240 ± 15 


275 ± 24 


7 


400 


55 ± 3 


117 ± 3 


144 ± 4 


179 ± 6 


206 ± 12 


241 ± 12 


7 


600 


55 ± 3 


101 ± 2 


123 ± 4 


156 ± 5 


183 ± 7 


222 ± 8 


6 


290 


48 ± 3 


161 ± 4 


190 ± 7 


223 ± 13 


272 ± 16 


314 ± 27 


6 


400 


55 ± 3 


144 ± 4 


170 ± 5 


222 ± 7 


258 ± 14 


282 ± 13 


6 


600 


56 ± 3 


124 ± 3 


152 ± 5 


200 ± 6 


244 ± 9 


271 ± 9 


5 


290 


11 ± 4 


75 ± 5 


100 ± 4 


137 ± 8 


151 ± 10 


184 ± 18 


5 


400 


22 ± 2 


79 ± 4 


98 ± 5 


126 ± 6 


143 ± 15 


161 ± 8 


5 


600 


26 ± 2 


65 ± 2 


81 ± 3 


108 ± 4 


123 ± 6 


153 ± 7 


<5 


290 


33 ± 16 


414 ± 21 


692 ± 45 


1100 ± 99 


1660 ± 142 


2239 ± 260 


<5 


400 


51 ± 15 


441 ± 21 


728 ± 26 


1225 ± 47 


1700 ± 99 


2317 ± 129 


<5 


600 


52 ± 14 


475 ± 13 


708 ± 36 


1205 ± 48 


1727 ± 66 


2270 ± 62 



TABLE IX. Fragment production cross section for a 24 Mg beam on elemental targets. The beam energy was 400 MeV/nucleon 
at extraction. Cross sections are in millibarns. 



■^frag 


H 


C 


Al 


Cu 


Sn 


Pb 


ll 


72 ± 3 


122 ± 3 


154 ± 5 


195 ± 8 


232 ± 14 


274 ± 22 


10 


63 ± 3 


112 ± 3 


135 ± 4 


176 ± 8 


197 ± 12 


250 ± 19 


9 


27 ± 1 


56 ± 2 


64 ± 3 


86 ± 4 


96 ± 7 


96 ± 9 


8 


56 ± 3 


113 ± 3 


136 ± 4 


175 ± 7 


181 ± 11 


232 ± 18 


7 


34 ± 2 


85 ± 2 


114 ± 4 


142 ± 6 


180 ± 11 


180 ± 14 


6 


33 ± 2 


113 ± 3 


150 ± 5 


200 ± 8 


191 ± 12 


280 ± 21 


5 


12 ± 3 


53 ± 4 


71 ± 5 


97 ± 8 


120 ± 12 


132 ± 16 


< 5 


31 ± 18 


374 ± 18 


656 ± 38 


1173 ± 79 


1597 ± 79 


2283 ± 239 



going from 290 to 400 and 600 MeV/nucleon, with ra- 
tios of 2.09 ± 0.09, 2.22 ± 0.08, and 2.32 ± 0.06 for the 
three energies, respectively. The results for 20 Ne show 
the same trend, although all ratios are smaller for the 
lighter beam. (Data points with 20 Ne at higher energies 
would be of interest.) For 24 Mg at 400 MeV/nucleon, we 
find a ratio of 2.32 ± 0.05, compatible with that found 
for 28 Si at the same beam energy. 



VII. SMALL ACCEPTANCE SPECTRA: 
GENERAL CONSIDERATIONS 

In the following, we again present the results grouped 
by beam ion, similar to the way in which the large- 
acceptance results were presented above. However, here 
it will be seen that the differences between corrected re- 
sults at different beam energies are, in many instances, 
quite large, for reasons that are related to the variations 
in the acceptances (which depend on beam ion, energy, 
and experimental configuration) and the associated cor- 
rection factors. Before presenting the results, discussion 



of some broadly-applicable points is in order, as there are 
several caveats in the interpretation of these data. The 
corrected light fragment cross sections presented here 
must be interpreted with considerable caution. 



A. Events with No Detected Fragments 

In some runs, the most populated region of the small- 
acceptance spectrum is that below the peak due to sin- 
gle charge-1 particles. These are events in which either 
nothing at all was detected, or the detected particle(s) 
deposited less than a minimum-ionizing charge 1 par- 
ticle. These events may contain low-energy electrons 
(£ w 1 MeV), or Compton electrons produced in the 
detector from the traversal of an energetic gamma-ray. 
These events are counted in the sense that they con- 
tribute to the charge-changing cross section, as invari- 
ably some charge is recorded at large acceptance, but no 
corrections are applied to this event category. We do not 
report these cross sections in the following. 
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B. Non-leading Charge-1 Fragments 

Interpretation of the cross sections for the lightest 
species, H and He, must be approached with particu- 
lar caution. These fragments are copiously produced in 
the interactions studied here, but typically they are not 
likely to be the highest- Z fragment detected in a given 
event, even at small acceptance. The experimental ap- 
proach used here cannot fully account for all of these 
fragments and therefore the cross sections obtained be- 
fore acceptance corrections are applied represent the de- 
tected portion of the overall production cross sections. 
Given the experimental setup, it is only feasible to mea- 
sure the cross section for events in which either a H or 
He ion is the highest-charge fragment seen in the small- 
acceptance detector, and extrapolate these results using 
our angular acceptance model. Even with these correc- 
tions, we expect that large portions of the true produc- 
tion cross sections are missed. Consider the simplest case 
in which a single proton or deuteron is stripped from the 
projectile (i.e., AZ — 1), a hydrogen fragment must also 
be created; however, the probability that it is detected 
is small, for two reasons. First, the heavier fragment an- 
gular distribution is much more sharply forward-peaked 
than that of the H fragment. Conservation of transverse 
momentum dictates that even in those rare events where 
the heavier fragment is produced at a large enough angle 
to miss the small-acceptance detector, the H fragment or 
fragments (for AZ > 1) must have balancing transverse 
momenta, so that they tend to be far outside the small 
acceptance. Secondly, the heavier fragment will virtually 
always be seen in the large-acceptance detector and has a 
very high probability to be seen in the small-acceptance 
detector, and with unsegmented detectors having limited 
resolution in AE, a non-leading H fragment will typically 
not be detected even if it is within the acceptance be- 
cause its contribution to the effective detected charge is 
so small. For example, in oxygen beam data, we cannot 
distinguish between the detection of a nitrogen fragment 
alone, Z^ s = 49, and coincidental detection of a nitrogen 
fragment and a proton, Z% s = 50. (By Z e g, we mean the 
"effective" Z seen in a particular detector pair; as a prac- 
tical matter, this is proportional to the square root of AE 
and - for ions with equal velocities - also proportional to 
the square root of the sum of the charges squared. See 
Ref. [7[ for additional explanation.) 



C. Non-Leading He Fragments 
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FIG. 4. Part of the charge spectrum from 1 GeV/nucleon 16 
ions incident on a polyethylene target of 2.82 g cm -2 depth. 
The H and He fragments have, on average, higher velocities 
than the primaries used to set the charge scale, and so these 
peaks appear below the corresponding integer values. Other 
peaks, or shoulders on larger peaks, are due to non-leading 
fragments and appear at effective charges of ~ 1.3 (two H 
fragments), 2.1 (He + H), etc. 



tra, as was shown in Figure 2 for a 14 N beam and again, 
with some important differences, in Figure 4 for a 16 O 
beam. As discussed previously, the peak around charge 
3.45 is due to the detection of three He fragments in co- 
incidence (although a contribution from Li + He cannot 
be ruled out). The peak in the charge 2.6 to 3.2 region 
is due to events in which pairs of He fragments are de- 
tected in coincidence, along with events in which there 
is a leading Li fragment. In contrast to the spectrum in 
Figure 2, here we do not see a clear separation of the Z 
pa 2.8 peak from the Li peak, and there also appears to 
be a comparatively larger share of events near 2.8. As 
in Figure 2, the peak near 1.3 is due to the detection of 
two charge-1 fragments in coincidence, and here a small 
peak is seen near charge 1.6 (in the low-end tail of the 
He distribution), which is likely due to detection of three 
charge-1 fragments in coincidence. Finally, the peak seen 
near effective charge of 2.1 is almost certainly due to co- 
incidences of helium and charge-1 fragments (as is the 
similar peak in Figure 2). 



Many of the same points made in the preceding sub- 
section apply to non-leading He fragments, but with a 
few important differences: (1) the He fragment angular 
distributions are more forward-focused than those of H 
fragments, and (2) they make larger contributions to ef- 
fective charge than do H fragments. Accordingly, the 
presence of non- leading and/or multiple He fragments 
can often be inferred from typical small-acceptance spec- 



We do not count non-leading fragments as contribut- 
ing to the production cross sections of a given species. To 
properly perform that analysis would require a more so- 
phisticated experiment. However, we will partially take 
account of the non-integer peaks in the following by mak- 
ing use of the events in the peaks due to coincidences of 
two or three helium fragments. 
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D. Lithium/Helium-Pair Cross Sections 

In the region around Zf lag = 3, there are contributions 
from events with either a leading lithium fragment or a 
coincidence of two helium fragments. Cross sections have 
been obtained by applying the average of corrections for 
fragment masses of 6 and 7 amu. Again, there is ambi- 
guity here; the "charge 3" peaks in all likelihood contain 
significant contributions from 8 Be (leading to a pair of 
4 He fragments), 7 Li, 6 Li, and pairs of 4 He fragments that 
are produced independently. Choosing to average A = 6 
and A = 7 represents a best guess as to the midpoint of 
the acceptance for these events. 



E. Three-Helium Fragment Production Cross 
Sections 

Although acceptance corrections are applied in all 
other cases, we choose not to apply them to the cross 
sections corresponding to the Z « 3.5 peaks. This is be- 
cause we have no a priori knowledge of the angular dis- 
tributions of the fragments observed in these events. It is 
possible to treat the three fragments as having been pro- 
duced independently, i.e., with no mutual correlations, in 
which case the corrections tend to be large. However, it 
may be that a significant fraction of these events arises 
from production of 8 Be in conjunction with another he- 
lium fragment; the decay of 8 Be produces two 4 He frag- 
ments that together have a more forward-focused angular 
distribution than a single He fragment. Because there is 
negligible Q in the 8 Be decay, the initial forward-going 
trajectory of the fragment tends to be preserved by the 
two 4 He decay products. The probability for detection 
of all three He fragments also depends on the Coulomb 
multiple scattering the fragments undergo, which occurs 
independently. The detection efficiency for these events is 
therefore a strong function of the production mechanism, 
the beam ion and energy, target material and depth, and 
the acceptance angle of the small-acceptance detectors 
used in the analysis. Of course, cross sections are by 
definition independent of the detection efficiencies. How- 
ever in this instance, uncertainties about the production 
mechanism make it impossible to estimate angular accep- 
tances with confidence, and we therefore present these 
cross sections without acceptance corrections. We will 
return to this subject repeatedly in the following as re- 
sults from each beam ion are presented and discussed. In 
the subsequent discussion, it will be convenient to define 
two hypotheses: (1) the three 4 Hc fragments produced 
independently, so that the detection efficiency, e, is sim- 
ply the cube of the efficiency for detection of a single 4 He 
fragment, i.e., £3^f e =[e( 4 He)] 3 ; (2) two of the three he- 
lium fragments are the products of 8 Be decay, and the 
third is produced independently, so that £3ij e =[e( 4 He)] 
x [e( 8 Be)]. We will refer to these as hypothesis 1 and 
hypothesis 2. 



F. Beryllium Production Cross Sections 

The charge 4 category also presents some ambiguities. 
For the 14 N beam, there is no plausible background from 
the detection of four helium fragments in coincidence 
(Z c g — 16). However, for all other beams reported on 
here, such events are at least theoretically possible, and 
if any such events occur they are indistinguishable from 
events in which a single Be fragment is detected. For 
the 16 beam ions, fission into a pair of 8 Be fragments, 
and subsequent decay of the 8 Be seems to be a plau- 
sible source of background. In the other direction, the 
detected Be cross section cannot include any direct con- 
tribution from the production of a single 8 Be fragment. 
When comparing to model calculations, it is necessary to 
subtract the predicted 8 Be contribution to the total Be 
production cross section. 

VIII. LIGHT FRAGMENT CROSS SECTIONS 

In the following, we present the cross sections obtained 
with the "small" acceptances shown in Table I above. 
Since the acceptances vary with beam ion, energy, and 
the angle subtended by the small acceptance detector in 
a given experiment, in each subsection we show a plot 
of the detection efficiency vs. fragment mass number for 
each of the experiments discussed therein. When apply- 
ing corrections, we assume charge 4 fragments have mass 
8, charge 2 fragments have mass 4, and charge 1 frag- 
ments have mass 1. For charge 3, we average the results 
obtained for masses 6 and 7. The choice of mass 8 for 
charge 4 represents a rough average of the acceptances for 
the stable isotopes with masses 7 and 9. The results are 
somewhat sensitive to these choices, but in the absence 
of isotopic resolution in the experiments, we are guided 
by the NUCFRG2 model, since (unlike LAQGSM and 
PHITS) the code directly outputs isotopic cross sections. 

A. Nitrogen Beams 

Table X shows the results for the production of Be and 
lighter fragments with the 14 N beams. The results for 
Be fragments were shown in Table VI above but are re- 
peated here as they help to illuminate the overall trends. 
As shown in Table I, the 290 MeV/nucleon data were 
obtained with an acceptance angle of 1.7° and the 400 
MeV/nucleon data at 3.8°. Figure 5 shows the calcu- 
lated acceptances as a function of fragment mass num- 
ber for the two nitrogen beams in the small-acceptance 
detectors. The larger efficiency in the 400 MeV/nucleon 
experiment is due both to the higher beam energy and to 
the larger acceptance angle. The acceptance corrections, 
applied to both data sets (except for the three-helium 
coincidence results) bring the charge 4 results into rea- 
sonable agreement at the two energies; the small differ- 
ences seen for charge 4 could well be real, that is, due 
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to actual weak energy dependence of the cross sections, 
although the differences are for the most part within or 
barely beyond one-sigma significance. 

The uncorrected results for coincidences of three he- 
lium fragments are, in contrast, very different from one 
another - the 400 MeV/nucleon cross sections are larger 
by a factor of about 3 for all targets. This is due to the 
combined effects of the larger acceptance angle and the 
greater forward focusing of fragments produced by the 
higher-energy beam. As described above, there is ambi- 
guity about the production mechanism for these events, 
and it is not clear how the cross sections should be cor- 
rected. We find, suggestively, that in our acceptance 
model for three helium fragments, assuming no energy- 
dependence of the cross sections, the ratio of efficiencies 
predicted by hypothesis (1) for these beams is a little 
over 5, and by hypothesis (2), exactly 3.0 [28| . 

For the Li/He-pair category, the results for the two 
energies are mutually consistent within the uncertain- 
ties. There is no particular reason to expect this. Con- 
sider that the two-helium contribution to the peak likely 
consists of events in which there were actually three He 
fragments produced, but only two were detected, so there 
is "feed-down" from the three-helium category. There is 
also feed-down to the next category (one leading He frag- 
ment detected) - events in which two He fragments are 
produced, but only one is detected. Given this complex- 
ity, and the differences in acceptances in the two exper- 
iments, it is surprising (and perhaps not meaningful) to 
find agreement. 

For the single-helium category, the acceptance- 
corrected cross sections at 290 MeV/nucleon are all 
larger than those at 400 MeV/nucleon. The ratios (400 
MeV/nucleon cross sections divided by 290 MeV/nucleon 
cross sections) are all mutually consistent, in the range 



0.70 - 0.80, for carbon and heavier targets. The fact 
that cross sections in this category appear larger at 290 
MeV/nucleon is not surprising. At the lower energy, 
there is more feed-down from the two- and three-helium 
categories due to more fragments being outside the ac- 
ceptance, and the acceptance corrections are substan- 
tially larger. Similar arguments apply to the Z = 1 
category, where again the cross sections obtained at 290 
MeV/nucleon are larger than those at 400 MeV/nucleon. 
For carbon and heavier targets, ratios are again mutually 
consistent, in the 0.78 to 0.93 range. 



B. Oxygen Beams 

Table XI shows the results for Be and lighter frag- 
ments with 16 beams at four energies. The corrected 
Be cross sections show no consistent pattern of energy 
dependence. For Sn targets, there is an increase with 
increasing beam energy, but this is not seen for other 
targets. 

To understand the cross sections in the other cate- 
gories, it is useful to note that for the 290, 600, and 1000 
MeV/nucleon experiments, the small-acceptance detec- 
tor angles were comparable, 1.7°, 1.4°, and 1.4°, respec- 
tively, while for the 400 MeV/nucleon experiment, the 
small acceptance detector subtended a half-angle of 2.5°. 
The calculated acceptances are shown in Figure 6. The 
curves for the 290 and 600 MeV/nucleon experiments sit 
almost on top of one another; this is fortuitous, and not 
by design. 

The effects of the acceptance differences can be seen in 
the uncorrected cross sections for three helium fragments. 
Results with the 290 and 600 MeV/nucleon beams are 
mutually consistent, which makes sense in view of the 
nearly-identical acceptance curves for the two experi- 
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TABLE X. Light fragment production cross section for N beams on elemental targets. Beam energies at extraction are shown, 
in units of MeV/nucleon. All cross sections are in millibarns. The results for three helium fragments detected in coincidence 
are not corrected for acceptance. 



■^frag 


■S'bcam 


H 


C 


Al 


Cu 


Sn 


Pb 


4 


290 


14 ± 2 


48 ± 2 


67 ± 3 


97 ± 5 


115 ± 9 


133 ± 11 


4 


400 


15 ± 2 


43 ± 2 


59 ± 4 


81 ± 5 


92 ± 7 


130 ± 9 


3 He coin. 


290 


8 ± 1 


18 ± 1 


25 ± 2 


28 ± 2 


31 ± 4 


28 ± 4 


3 He coin. 


400 


27 ± 3 


54 ± 3 


69 ± 4 


80 ± 5 


91 ± 7 


93 ± 7 


3 or 2 He coin. 


290 


49 ± 11 


138 ± 9 


186 ± 13 


296 ± 25 


319 ± 33 


436 ± 46 


3 or 2 He coin. 


400 


48 ± 7 


147 ± 7 


177 ± 10 


268 ± 14 


330 ± 20 


397 ± 23 


2 


290 


97 ± 26 


337 ± 21 


482 ± 32 


776 ± 63 


841 ± 83 


1222 ± 123 


2 


400 


53 ± 11 


267 ± 13 


345 ± 19 


543 ± 26 


673 ± 37 


927 ± 48 


1 


290 


92 ± 29 


399 ± 26 


712 ± 48 


1194± 98 


1549 ± 153 


2330 ± 232 


1 


400 


26 ± 15 


401 ± 19 


569 ± 31 


915± 45 


1211 ± 66 


1839 ± 81 



TABLE XI. Light fragment production cross section for O beams on elemental targets. Beam energies at extraction are 
shown, in units of MeV/nucleon. All cross sections are in millibarns. The results for three helium fragments detected in 
coincidence are not corrected for acceptance. 



■^frag 


■^bcam 


H 


C 


Al 


Cu 


Sn 


Pb 


4 


290 


8 ± 2 


49 ± 3 


66 ± 4 


71 ± 5 


68 ± 8 


114 ± 16 


4 


400 


16 ± 3 


46 ± 4 


61 ± 4 


75 ± 6 


87 ± 8 


98 ± 9 


4 


600 


14 ± 3 


45 ± 4 


53 ± 3 


66 ± 5 


96 ± 8 


112 ± 11 


4 


1000 


16 ± 3 


50 ± 4 


54 ± 3 


76 ± 4 


109 ± 7 


117 ± 9 


3 He coin. 


290 


5 ± 1 


17 ± 1 


20 ± 1 


17 ± 2 


26 ± 3 


22 ± 4 


3 He coin. 


400 


15 ± 3 


47 ± 4 


56 ± 4 


61 ± 5 


64 ± 6 


55 ± 5 


3 He coin. 


600 


7 ± 1 


20 ± 2 


20 ± 1 


23 ± 2 


26 ± 3 


29 ± 3 


3 He coin. 


1000 


15 ± 2 


33 ± 2 


41 ± 2 


48 ± 3 


48 ± 3 


59 ± 5 


3 or 2 He coin. 


290 


30 ± 6 


125 ± 7 


143 ± 11 


193 ± 15 


210 ± 20 


271 ± 28 


3 or 2 He coin. 


400 


36 ± 9 


136 ± 10 


168 ± 10 


205 ± 16 


240 ± 19 


308 ± 23 


3 or 2 He coin. 


600 


33 ± 6 


107 ± 8 


118 ± 6 


161 ± 9 


209 ± 13 


225 ± 17 


3 or 2 He coin. 


1000 


46 ± 7 


129 ± 9 


160 ± 8 


220 ± 11 


267 ± 15 


294 ± 19 


2 


290 


67 ± 15 


357 ± 18 


474 ± 34 


608 ± 43 


702 ± 60 


904 ± 92 


2 


400 


45 ± 15 


270 ± 19 


335 ± 20 


461 ± 33 


582 ± 43 


820 ± 56 


2 


600 


103 ± 20 


410 ± 29 


497 ± 24 


592 ± 31 


961 ± 54 


1224 ± 78 


2 


1000 


70 ± 14 


278 ± 19 


377 ± 17 


501 ± 25 


695 ± 37 


819 ± 49 


1 


290 


34 ± 16 


385 ± 20 


607 ± 44 


781± 57 


1033 ± 90 


1426 ± 148 


1 


400 


21 ± 17 


325 ± 23 


525 ± 31 


806± 57 


970 ± 71 


1471 ± 101 


1 


600 


97 ± 25 


515 ± 38 


719 ± 35 


967± 52 


1487 ± 86 


2017 ± 128 


1 


1000 


54 ± 18 


404 ± 28 


576 ± 27 


836± 42 


1274 ± 68 


1585 ± 94 



ments, with the additional assumption that the produc- 
tion cross section is weakly or not at all dependent on 
beam energy. For the 600 and 1000 MeV/nucleon data, 
where the acceptance angle was the same, the higher- 
energy beam yields larger cross sections in all cases, typ- 
ically by a factor close to 2, due to the greater forward 
boost, which again is reflected in the acceptance curve in 
Figure 6. Compared to results at other energies, the 400 
MeV/nucleon cross sections are significantly larger than 
those obtained at 290 and 600 MeV/nucleon, by roughly 
a factor of 2, and in most cases are larger than those 
obtained at 1000 MeV/nucleon. This can only be due 
to the comparatively large acceptance angle employed in 
the 400 MeV/nucleon experiment; the curves in Figure 
6 bear this out. The acceptance model predicts that, 
for hypothesis 1, the cross section measured in the 290 



MeV/nucleon experiment should be only 25% as large as 
that measured in the 400 MeV/nucleon experiment. For 
hypothesis 2, 8 Bc + 4 He, the same ratio is predicted to 
be 43%, somewhat closer to (but slightly higher than) the 
observed ratios. These ratios suggest that the observed 
events are a mix of the two types, with a majority of 
8 Be events. Comparing the 600 and 1000 MeV/nucleon 
results to those obtained at 400 MeV/nucleon yields a 
muddled picture: the 600 MeV/nucleon data are also 
consistent with a mix of the event types, but the 1000 
MeV/nucleon data are consistent with being entirely due 
to 8 Be + 4 He. 

Turning to the Li or 2 He cross sections, we find a 
high degree of consistency between the 400 and 1000 
MeV/nucleon experiments. At first glance, this is slightly 
surprising since the acceptances in the 400 MeV/nucleon 
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experiment are larger by 10-20%, as can be seen in 
Figure 6. The cross sections obtained with the 290 
MeV/nucleon beam are marginally consistent with, but 
in all cases smaller than, those obtained at 400 and 1000 
MeV/nucleon, and the 600 MeV/nucleon cross sections 
are (for the most part) smaller still. 

Based on the acceptance curves, we might expect the 
600 MeV/nucleon cross sections to be consistent in all 
cases with those from the 290 MeV/nucleon beam, but 
this is not what is observed. For the single-He category, 
the cross sections obtained at 600 MeV/nucleon are gen- 
erally the largest. Figure 7 shows the cross sections as 
functions of target mass for the (leading) Z = 2 and Z = 1 
categories. For both, the 600 MeV/nucleon cross sections 
are larger than those at the other energies, particularly 
for the Sn and Pb targets. In view of the similarity in 
the detection efficiencies between this experiment and the 
290 MeV/nucleon experiment, this may seem surprising, 
but referring to Tables II and VII is illuminating. The 
charge-changing cross sections tend to be slightly larger 
at 600 MeV/nucleon compared to 290 MeV/nucleon (Ta- 
ble II), and the cross sections in the Z < 5 category are 
10-20% larger at the higher energy (Table VII). Small 
differences in the production cross sections are multiplied 
by large acceptance corrections to produce the large dif- 
ferences seen in Figure 7. And although the grouped 
Z< 5 production cross sections are found to be approxi- 
mately equal at 600 and 1000 MeV/nucleon (Table VII), 
the higher-energy results receive significantly smaller cor- 
rections. 



compared to 1.7°, which is as expected. Taking the ratio 
of the cross sections at 2.5° to those at 1.7° by target, 
we see values ranging from 1.2 to 2.2, but on closer in- 
spection we find the data are reasonably consistent with 
being independent of the target. The weighted average 
of the ratios is 1.57 ± 0.09, with a \ 2 of 7.3 for 5 degrees 
of freedom for the ratio being target-independent. Our 
acceptance model predicts, using the same data points 
that populate Figure 8, a ratio of 1.74 for hypothesis 2 
and a ratio of 2.82 for hypothesis 1. The data are clearly 
- again - much more consistent with hypothesis 2 ( 8 Be + 
4 He). At 290 and 400 MeV/nucleon, the cross sections 
in this category are (except for the hydrogen target) sig- 
nificantly larger at the higher energy. This is due to the 
larger acceptance in the 400 MeV/nucleon experiment 
illustrated in Figure 8. 

For the Li or He pair category, there is a consistent 
ordering of the cross sections, with those obtained at 290 
MeV/nucleon being the largest for all targets and those 
obtained at 600 MeV/nucleon being the smallest. The 
400 MeV/nucleon data are in between, but are mostly 
close to the 600 MeV/nucleon results. Referring to Ta- 
ble VIII, the (uncorrected) cross sections for the grouped 
Z< 5 category are seen to be mutually consistent within 
uncertainties, suggesting that the ordering of the cross 
sections in this category is likely due to the acceptance 
corrections and not to any energy dependence in the un- 
derlying production mechanism. Similar remarks apply 
to the Zf rag = 2 and 1 categories, although in the lat- 
ter we are unable to report cross sections from the 290 
MeV/nucleon experiment. 



C. Neon Beams 

In Table XII, we show results for charge 4 and lighter 
fragments with 20 Ne beams at three energies. Accep- 
tance curves for the three beam energies are shown, for a 
2 g cm -2 12 C target, in Figure 8. There were two experi- 
ments with the 600 MeV/nucleon beam, one of which was 
reported on earlier Q. Small-acceptance cross sections 
shown here are the result of combining the data sets as 
described above, except for the three helium-fragment co- 
incidence category, for which we show separately results 
obtained with acceptances of 2.5° and 1.7°. 

The charge 4 production cross sections increase with 
increasing energy for the hydrogen target, but for other 
targets tend to decrease slightly with increasing energy. 
This is quite similar to the behavior seen for smaller 
charge changes in Tables VII and VIII. The cross sec- 
tions for three helium fragments detected in coincidence 
all increase with energy, but this may be largely an ac- 
ceptance effect. Since we do not know a priori whether 
or not there is any dependence of the underlying pro- 
duction cross section on beam energy, the fact that we 
have measurements at two acceptance angles using the 
600 MeV/nucleon beam is potentially instructive. For 
all targets, the cross sections for 600 MeV/nucleon beams 
are found to be larger when the acceptance angle is 2.5° 



D. 400 MeV/nucleon 24 Mg 

Table XIII shows the cross sections obtained for the 
400 MeV/nucleon 24 Mg beam, and Figure 9 shows the 
results of the acceptance calculation for this beam on 
a 4 g cm" 2 12 C target as was used in the experiment. 
Also shown in Figure 9 are the efficiency curves found 
for the other experiments with 400 MeV/nucleon beams 
reported here. The inset figure zooms in on the frag- 
ment charge range from 1 to 4, which we examine in the 
following. Again, the acceptances are only weakly de- 
pendent on the type and depth of the target material, so 
the curves in Figure 9 are representative of all runs with 
these beams. Note that if all other variables are held 
constant except for the mass of the beam ion, as that 
mass increases the acceptance of the lightest fragments 
decreases, because the width of the angular distribution 
is driven by a term proportional to (Abeam — ^frag) 1 ^ 2 - 

Since data are only available with a single beam energy 
for 24 Mg, we cannot probe these results for possible en- 
ergy dependence of the production cross sections. How- 
ever, other features of the data merit comment. For the 
hydrogen target, the production cross sections for both 
Zfrag = 1 and Zf rag = 2 are consistent with 0. This is 
seen for Zf rag = 1 for several of the other beams reported 
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FIG. 7. Corrected cross sections for events with a single leading charge 2 fragment (left) and a leading charge 1 fragment (right) 
as measured at small acceptance detectors for the four le O experiments. 



TABLE XII. Light fragment production cross section for Ne beams on elemental targets. Beam energies at extraction are 
shown, in units of MeV/nucleon. All cross sections are in millibarns. The results for three helium fragments detected in 
coincidence are not corrected for acceptance. For the 290 MeV/nucleon beam, fragments with charge 1 could not be resolved. 



■^frag 


-E'bcam 


H 


C 


Al 


Cu 


Sn 


Pb 


4 


290 


6 ± 1 


52 ± 4 


66 ± 5 


96 ± 8 


130 ± 14 


135 ± 17 


4 


400 


9 ± 2 


46 ± 3 


59 ± 3 


70 ± 4 


82 ± 7 


103 ± 9 


4 


600 


14 ± 2 


46 ± 3 


59 ± 4 


68 ± 5 


102 ± 7 


96 ± 7 


3 He coin. 


290 


8 ± 1 


14 ± 2 


21 ± 1 


13 ± 2 


16 ± 3 


14 ± 4 


3 He coin. 


400 


8 ± 1 


21 ± 1 


25 ± 1 


31 ± 1 


35 ± 3 


30 ± 3 


3 He coin. 


600 (2.5°) 


16 ± 2 


56 ± 1 


58 ± 2 


70 ± 5 


112 ± 6 


81 ± 4 


3 He coin. 


600 (1.7°) 


14 ± 3 


31 ± 1 


37 ± 2 


53 ± 3 


51 ± 5 


55 ± 5 


3 or 2 He coin. 


290 


27 ± 2 


161 ± 10 


257 ± 17 


371 ± 28 


294 ± 28 


632 ± 66 


3 or 2 He coin. 


400 


20 ± 6 


119 ± 6 


149 ± 7 


191 ± 10 


229 ± 15 


272 ± 18 


3 or 2 He coin. 


600 


19 ± 7 


103 ± 8 


128 ± 10 


163 ± 14 


219 ± 19 


262 ± 22 


2 


290 


73 ± 4 


361 ± 22 


526 ± 33 


973 ± 65 


1293 ± 106 


1632 ± 150 


2 


400 


40 ± 14 


298 ± 15 


402 ± 19 


547 ± 28 


604 ± 37 


881 ± 53 


2 


600 


18 ± 11 


190 ± 15 


265 ± 21 


378 ± 32 


434 ± 37 


588 ± 47 


1 


400 


8 ± 18 


392 ± 19 


572 ± 26 


809± 51 


998 ± 65 


1476 ± 96 


1 


600 


11 ± 14 


247 ± 19 


455 ± 36 


680± 57 


871 ± 72 


1232 ± 99 



above, but not for Zf lag — 2. The large uncertainties 
in these data do not permit us to make a stronger state- 
ment; comparing to the 400 MeV/nucleon 20 Ne result for 
Zfmg = 2, we can only say that the cross section with the 
24 Mg beam is smaller at about the 1.5a level. It would 
be useful to obtain more data with this beam ion, at 400 
MeV/nucleon and other energies, particularly since Mg 
is prominent among the GCR heavy ions. 

Comparing these data to the 400 MeV/nucleon 20 Ne 
cross sections, we find very similar values for the Zf ra g 
— 4 cross sections for all targets, and also for the three- 
helium-coincidence category. For the lower- .Z categories, 
excluding the hydrogen target data, the cross sections 
for the 24 Mg beam are significantly larger than for the 
20 Ne beam, and also for the other beam ions at 400 



MeV/nucleon. The same can also be said for the Zf lag 
— 1 category. Figures 10 and 11 illustrate these trends. 
Cross sections for the beams other than Mg are in many 
instances mutually consistent, with the 14 N beam data 
tending to be slightly higher than the 16 and 20 Ne data 
for Zf r ag = 1. The 24 Mg cross sections are larger than 
those for other beams by factors of about 1.5 for Zf rag — 2 
and factors of about 2 for Zf lag = 1 . Referring to the inset 
figure in Fig. 9, we see that the detection efficiencies for 
these light fragments were lowest in the Mg experiment. 
We expect that production cross sections for protons and 
helium must increase with increasing beam charge and 
mass since there are simply more nucleons available to 
be sheared off the projectile, and the probability for de- 
tecting one or more fragments should increase. 
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TABLE XIII. Light fragment production cross section for a 400 MeV/nucleon 24 Mg beam on elemental targets. All cross 
sections are in millibarns. The results for three helium fragments detected in coincidence are not corrected for acceptance. 



■^frag 


H 


C 


Al 


Cu 


Sn 


Pb 


4 


7 ± 2 


41 ± 3 


55 ± 4 


88 ± 8 


78 ± 9 


126 ± 16 


3 He coin. 


6 ± 1 


15 ± 1 


25 ± 2 


28 ± 3 


28 ± 3 


31 ± 5 


3 or 2 He coin. 


17 ± 7 


139 ± 9 


183 ± 13 


264 ± 20 


274 ± 25 


321 ± 35 


2 


± 19 


428 ± 27 


586 ± 40 


783 ± 57 


1074 ± 89 


1470 ± 138 


1 


25 ± 27 


568 ± 37 


1109 ± 75 


1957 ± 140 


2545 ± 209 


4264 ± 389 
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FIG. 8. Acceptance vs. fragment mass number as calculated 
for the small acceptance detectors for the four 20 Ne experi- 
ments. 
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FIG. 10. Acceptance-corrected production cross sections for 
the leading helium category for the 400 MeV/nucleon beams. 
Hydrogen-target data are excluded. 
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FIG. 9. Acceptance vs. fragment mass number as calculated 
for the small acceptance detectors for the experiments with 
400 MeV/nucleon beam energies and carbon targets. 
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FIG. 11. Acceptance-corrected production cross sections for 
the leading Z = 1 category for the 400 MeV/nucleon beams. 
Hydrogen-target data are excluded. 



IX. FRAGMENT CROSS SECTION 
COMPARISONS 



Articles in the literature for the beams studied here 
(and similar beams) have reported fragment cross sec- 
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tions for small charge changes, limiting the scope of 
model comparisons. Here, we are able to extend the com- 
parisons down to Z = 2, albeit with somewhat large un- 
certainties owing to the nature of the measurements. We 
do not believe the Z — 1 data reported above are suit- 
able for model comparisons because there is no feasible 
way to perform multiplicity weighting, mostly because a 
dominant share of the charge 1 fragments produced are 
undetectable since they accompany much heavier frag- 
ments. 



A. Model Descriptions of Fragmentation Processes 

Of the three models used here, PHITS and LAQGSM 
are three-dimensional Monte Carlo codes that simu- 
late individual interactions in great detail. In contrast, 
NUCFRG2 is an engineering code designed for computa- 
tional speed; its physics content is based on an abrasion- 
ablation formulation with its free parameters tuned to 
reproduce p-nucleus cross section data. 

In LAQGSM, when the mass number of the excited nu- 
cleus is greater than 12, a three-stage process is modeled: 
intranuclear cascade (INC); preequilibrium emission of 
fragments from the excited remnant nucleus; and evapo- 
ration and/or fission, if the compound nucleus is heavy 
enough to fission. When the mass number of the excited 
nucleus is 12 or less, LAQGSM uses the Fermi break-up 
model after INC. That is, the Fermi breakup model is 
used before the preequilibrium stage, as well as during 
the preequilibrium and evaporation stages, if the mass 
number of the excited nucleus becomes less than 13 due 
to emission of some preequilibrium and/or equilibrium 
particles. 

For the beam ions studied here, a significant share of 
the reactions fall into the latter category and therefore 
test the Fermi break-up model. The latest version of 
the LAQGSM (03.03) incorporates an improved version 
of the Dubna Cascade Model [2{|. It uses a continu- 
ous nuclear density distribution and experimental cross 
sections at energies below 4.5 GeV/nucleon. It has pre- 
viously been shown [2^, [2(| that there is good agreement 
between LAQGSM version 03.03 calculations and data 
taken by our group with 28 Si beams at 290, 600, and 
1200 MeV/nucleon [8j for fragment production cross sec- 
tions. 

PHITS uses a cross section model developed by Tri- 
pathi to determine interaction probabilities, the Jaeri 
Quantum Molecular Dynamics (JQMD) model to de- 
scribe nucleus-nucleus collisions, and the Generalized 
Evaporation Model (GEM) to model fission and evap- 
oration processes. Several previous comparisons between 
our data and PHITS can be found in the literature 

@,[I3H[I!11. 



B. Production Cross Sections for Charges 2 and 3 

In the preceding, we reported cross sections in three 
categories that pertain to Li and He production: "3 He" , 
"Li/2 He", and ll Z — 2." Cross sections obtained in 
the latter two categories are corrected for acceptance as- 
suming fragment mass numbers of 7 and 4, respectively. 
Multiplicity-weighted estimates of the total He produc- 
tion cross sections can be made with two opposing as- 
sumptions, in the hope of bracketing the true cross sec- 
tions. For a given beam ion/energy /target combination, 
our lower-bound estimate is obtained by multiplying the 
"3 He" cross section by three and adding it to the U Z 
= 2" cross section. The upper-bound estimated is ob- 
tained by multiplying the "Li/2 He" cross section by two 
and adding this to the lower-bound estimate. These esti- 
mates do not include contributions from non-leading He 
fragments that are undetectable with our experimental 
methods, so it is conceivable that even our upper-bound 
estimates are too small. Also, there is no acceptance 
correction for the "3 He" category, which would tend to 
cause an error in the same direction. An error in the 
other direction stems from the implicit assumption that 
there is no lithium production whatsoever, which cannot 
be true. 

On average, the upper-bound estimates give better 
agreement with the models, particularly for PHITS and 
LAQGSM. That is, the agreement is not necessarily bet- 
ter with the upper bound estimate for each data point, 
but lower cumulative uncertainty (defined below) is ob- 
tained for all three models when the predictions are com- 
pared to the upper-bound estimates. 

For the 14 N and 16 beams, the estimates of total 
He production appear to be largely energy-independent. 
Both upper- and lower-bound values for a given target 
tend to be mutually consistent across beam energies. 
Only the upper-bound values for hydrogen targets show 
a statistically-significant slope (a slight rise) when plot- 
ted as a function of the beam energy. The left-hand plot 
in Figure 12 shows the dependences of these estimates 
for 16 beams on H, C, and Al targets. Similar results 
are obtained for Cu, Sn, and Pb targets, but only the 
results for the three lighter targets are shown in order 
to keep the plot readable. In view of the varying exper- 
imental acceptances (see Figure 8 above), this is quite 
an interesting result. However, the 20 Ne upper and lower 
bounds for total helium production do not quite show the 
same behavior. As can be seen in the plot on the right- 
hand side of Figure 12, the results with the 400 and 600 
MeV/nucleon beams tend to be mutually consistent, but 
the cross section estimates with the 290 MeV/nucleon 
beam are in all cases significantly larger. 

To estimate the Li cross sections, we simply take half 
the measured cross section for the "Li/ 2 He" category 
and assign a relative error of ± 33%. This allows for the 
extreme (and implausible) possibilities that either all of 
these events are Li, or that none of them is Li, at the 
three-sigma level. 
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FIG. 12. Energy dependence of helium production cross sections as estimated by methods described in the text, H, C, and Al 
targets. Both upper- and lower-bound estimates from the data are shown. Left, in Fig. 12(a), results for le O beams; right, 
(b), for 20 Ne beams. In Fig. 12(a), the lines shown are for linear fits. 



Model comparisons to the Z = 1 data are excluded 
from consideration here for reasons mentioned above. 
Cross sections for Z = 1 predicted by PHITS and 
LAQGSM are typically factors of 2-3 larger than mea- 
sured cross sections. For NUCFRG2, the discrepancies 
are typically factors of 10-20, with the model predictions 
again being larger than the measured. 



C. Model Predictions of Helium Production by O 
and 20 Ne Beams 

Before proceeding to the overall comparisons of frag- 
ment cross sections and models, here we focus on the 
estimated helium production cross sections compared to 
the models. Given the difficulties involved in obtaining 
the estimates, the relatively large uncertainties associ- 
ated with them, and the fact that we do not include 
non- leading He fragments in our measurement, it is rea- 
sonable to wonder whether there is any correspondence 
at all with the model predictions. 

The cross sections in Tables X through XII that involve 
helium fragments are somewhat difficult to interpret, 
largely because the varying acceptances in the different 
experiments lead to disparate results for the different cat- 
egories. However, the method of bounding the produc- 
tion cross sections described above appears to at least 
partially compensate for the acceptance effects, based on 
the results shown in Figure 12 and on the reasonable 
agreement found with PHITS and LAQGSM predictions 
that we describe here. 

In Figure 13, we plot the upper-bound estimates of the 
He production cross sections from the data against the 
multiplicity-weighted cross sections predicted by PHITS 



and LAQGSM for the two 14 N beams and four 16 beams 
(left), and separately for the three 20 Ne beams and one 
24 Mg beam (right). The experimental error bars are 
shown. For the 16 beams incident on H, C, Al, and 
Cu targets (measured cross sections below 1100 mb), the 
agreement with PHITS is good, while the measured cross 
sections for the Sn and Pb targets are systematically 
higher than predicted by PHITS by about one standard 
error. The comparisons shown here are to the upper- 
bound estimates from the data, and looking at only the 
16 data, one might conclude that this is what causes 
almost all of the data points to fall above the 45° line. 
However, the situation is reversed for the 20 Ne beams: 
almost all of the predicted cross sections are larger than 
the measured cross sections, for both models, and the 
LAQGSM predictions are generally closer to the data 
than are the PHITS predictions. 

The good agreement between the measured upper- 
bound values and PHITS predictions for 16 O beams is 
remarkable, but given that a similar level of agreement 
is not seen with 20 Ne and 24 Mg beams, we conclude that 
it is likely fortuitous. However, in general, it is worth 
noting that the data and the predictions of the two mod- 
els are not grossly in disagreement. The dotted lines in 
Figures 13(a) and 13(b) show factor-of-two errors; in the 
left-hand plot, all points fall between the lines, and in 
the right-hand plot, five H-target data points fall outside 
the lines while all other points are within. 



D. Cumulative Relative Uncertainty 

Employing the bounding method for helium produc- 
tion cross sections, and taking the upper-bound esti- 
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Helium Cross Section, Model (mb) 

FIG. 13. Upper-bound estimates for helium production cross sections plotted vs. PHITS and LAQGSM predictions, for C, Al, 
Sn, and Pb targets. Results for le O beams are shown in Fig. 13(a) on the left, and for 20 Ne and 24 Mg beams in Fig. 13(b) 
on the right. For all beams, open symbols are used for LAQGSM and filled symbols for PHITS; in the left-hand plot, square 
symbols are for 16 O beams, and in the right-hand plot, square symbols are for the 24 Mg beam. The dotted lines correspond to 
factor-of-two differences between the data and the models. 



mates, we have 408 fragment production cross sections 
available for model comparisons. A concise method for 
comparing data and model predictions has recently been 
described by Norman and Blattnig [19|. Here, we briefly 
recapitulate the method and apply it to the ten data 
sets presented above. Comparisons are performed for 
NUCFRG2, PHITS, and LAQGSM for fragment species 
from one charge unit below the primary down to helium. 
The three quantities needed to implement the method 
are defined as follows. 

D+ = M ( Xi ) - [E (xj)+e(xj)] 
E(xi) + e(xi) 

D _ = MjxJ-lEjxJ-ejXi)] 
E (x^ - e (xi) 

U (an) = MAX (\D+ ( Xl )\ , \D~ ( Xi )\) 

The D value for a given data point i depends on the 
cross section predicted by the model, M(xi), the cross 
section as measured in the experiment, E(xi), and the 
uncertainty associated with the measurement, e(xi). Val- 
ues of D were calculated for each of the 408 data points 
available. For each model, for helium, comparisons to 
the upper-bound values described above yield the lowest 
average U(xi) values. 

Figure 14 shows the cumulative probability distribu- 
tions of the relative uncertainties as defined above, calcu- 
lated for NUCFRG2, PHITS, and LAQGSM. This graph 
should be read from left (0.0) to right (1.0). All three 
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FIG. 14. Cumulative relative uncertainty as defined in the 
text, calculated for the NUCFRG2, PHITS, and LAQGSM 
models. 



models yield similar curves for about 25% of the data, 
at which point the NUCFRG2 curve rises slightly above 
the other two. The PHITS and LAQGSM curves lie 
practically on top of one another, below the NUCFRG2 
curve, from about 0.25 to about 0.70. There, the 
LAQGSM curve begins to rise, crossing the NUCFRG2 
curve around 0.76. The LAQGSM curve remains above 
the NUCFRG2 curve, which remains slightly above the 
PHITS curve. The ordering persists out to about 0.99, 
where all three curves converge at a U (#,) value of about 
1.6. Based on these curves, PHITS can be said to have 
the best overall agreement with the data. 
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Interpreting the curves in Figure 14 is somewhat sub- 
jective. On one hand, we can say that with any of the 
models under consideration, about 70% of the cross sec- 
tions have U(xi) values less than 0.5. More critically, 
we can say that only about 30% of the predicted cross 
sections are highly accurate (U{xi) values less than 0.25) 
and that values of U{xi) above 1.0 represent significant 
weaknesses of the codes in those cases where the exper- 
imental errors are modest. (Large experimental errors 
dominate the U{xi) in a few cases here.) Of course the 
definition of what constitutes a "good" value of U(xi) is 
arbitrary and it could be argued that any values above, 
say, 0.5 are problematic. Figure 14 shows that, at least 
for these beams, PHITS is slightly more accurate overall 
than the other two models tested for these beam ions and 
energies. 



E. Model Systematics 

The cumulative relative uncertainty curves shown in 
Figure 14 give an overall picture of how well the models 
reproduce the data, but there are details and systematic 
differences that cannot be conveyed in such a compact 
format. A highly detailed discussion of these comparisons 
is beyond the scope of this article, but we will describe 
some obvious features that are noteworthy. We begin 
with LAQGSM, which, as shown in Figure 14, deviates 
from the other two models considered here once we get 
to about 70% of the data. A closer look at the source of 
the discrepancies reveals that they are largely confined 
to two fragment categories, charge 3, and AZ =1. To 
illustrate this, in Figure 15 we show the U(xi) values for 
LAQGSM predictions of N production cross sections vs. 
target mass for the four 16 beams. Clearly, LAQGSM 
does not agree well with the data for the lighter targets, 
and does somewhat better (but still not well) as target 
mass increases. In all cases, U{xi) is above 0.5, indicat- 
ing a lack of agreement with the data. The situation is 
similar, albeit not quite as bad, for AZ = 1 with the 
other beam ions. With the exception of the 14 N data, 
these cases test the INC part of LAQGSM. For Li pro- 
duction, which in LAQGSM depends on the Fermi break- 
up model, half of the values of U(xi) are above 1.0, and 
several more are above 0.9. While the U (xi) values for Li 
production tend to be comparatively large for all models 
due to the large experimental uncertainty, the LAQGSM 
comparison yields particularly large values, because (un- 
like the other models) it predicts Li cross sections that 
are in all cases significantly larger than the cross sections 
estimated from the data. With the unavoidable large 
uncertainties assigned to these data points, the D~ val- 
ues (corresponding to the case where the estimated cross 
section is smallest) are quite large. 

With sixty data sets (ten beams and six target ma- 
terials), the cross sections for Li and for AZ = 1 ac- 
count for 120 points of comparison out of 408. The lack 
of agreement for these cases largely explains the rise of 
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FIG. 15. Values of relative uncertainty U(xi) obtained with 
the LAQGSM model for the reaction 16 O + A -> N + X for 
the four beams and six targets used in the measurements. 



the LAQGSM curve in Figure 14 when the data fraction 
reaches 0.7. 

For NUCFRG2, about 30% of the cross sections pre- 
dicted by this model yield a U(xi) value above 0.5, as can 
be seen in Figure 14. The NUCFRG2 curve shows an in- 
flection point at a data fraction of about 0.94, followed by 
a fairly steep rise, largely driven by fourteen data points 
that have U(xi) values greater than 0.90. Ten of the 
fourteen discrepant data points are for hydrogen targets; 
five of these ten are for helium production. Of the other 
four highly discrepant points, three are for F production 
(■Zfi-ag = 9) in the 24 Mg beam data. In all fourteen cases, 
NUCFRG2 predicts cross sections that are roughly fac- 
tors of 2 larger than the measured cross sections. 

For PHITS, only about 20% of the cross sections have 
U(xi) values greater than 0.5, and only 5 of the 408 cross 
sections yield values above 0.90. Most of the large U{xi) 
values are seen are for Li production, where the experi- 
mental uncertainty contributes most of the total. 

All models yield large values of U(xi) for the helium 
production cross section with 24 Mg on H. The upper- 
bound estimate is 52 ± 20 mb, and the three model 
predictions are remarkably consistent for this point. 
NUCFRG2 predicts a cross section of 176 mb, LAQGSM 
178 mb, and PHITS 189 mb. All are more than a factor 
of three larger than the estimate from the data. This 
is the most extreme example of a broader trend for H 
targets, in which all three models predict substantially 
larger He production cross sections than are estimated 
from the data. 

The mean values of U(xi) for the 408 cross sections 
considered here are 0.35, 0.40, and 0.43 for PHITS, 
NUCFRG2, and LAQGSM, respectively. The order of 
the results (PHITS smallest, LAQGSM largest) is pre- 
dictable given the curves in Figure 14. Other compar- 
isons between PHITS, LAQGSM, and data have pre- 
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viously been reported, including measurements of neu- 
tron spectra [33j , proton-induced reactions [34|, fragmen- 
tation of 28 Si [35|, and fragmentation of 12 C and 56 Fe 
[26j . (Note that results reported as being for the MC- 
NPX code use LAQGSM.) Overall, LAQGSM tends to 
be at least as accurate as PHITS, and in some cases it is 
considerably more accurate. The results of the analysis 
presented here pertain to light-ion beam species; further 
investigation using cross section data for heavier beams 
is needed. 



F. Comparison to Webber et al. Data 

The Cumulative Relative Uncertainty method can also 
be used to compare a subset of the large-acceptance frag- 
ment production cross sections to those obtained by Web- 
ber et al. (lij . In the equations for D + and D~ above, we 
take the M(xi) terms to be the cross sections reported 
in the earlier work, the E(xi) to be our cross sections 
as per Tables VI through IX, and the e(xi) to be the 
quadrature sums of the uncertainties reported by the two 
groups. There are 54 data points that can be reasonably 
compared. We find an average U(xi) of 0.154 and no val- 
ues greater than 0.533. Although this level of agreement 
suggests lurking systematic errors beyond those claimed, 
it is still far better than the agreement between the data 
and any of the models. 

Fragment production cross sections with isotopic reso- 
lution using oxygen beams at comparable energies have 
been reported by Leistenschneider et al. [3al (at sa 
600 MeV/nucleon) and Momota et al. Hf] (at 290 
MeV/nucleon). However, comparisons to those experi- 
ments are not straightforward owing to large differences 
in the detector configurations and the methods of ex- 
tracting fragment production cross sections. The Leis- 
tenschneider et al. data have been compared to a previ- 
ous version of LAQGSM 

X. CONCLUSIONS 

We have presented 60 charge-changing cross sections 
and over 500 fragment production cross sections for 
beams of 14 N, 16 0, 20 Ne, and 24 Mg at energies rang- 
ing from 290 to 1000 MeV/nucleon. Comparisons to the 
PHITS, NUCFRG2, and LAQGSM models have been 
made. PHITS and NUCFRG2 can be directly compared 
to the measured charge-changing cross sections, and show 
reasonable agreement with the data. A simple geomet- 
ric model with two free parameters (the nucleon radius 
and the overlap term) can be tuned to fit the bulk of the 
charge-changing cross section data slightly better than 
the more sophisticated models. 

Large-acceptance charge spectra are used to extract 
fragment production cross sections for charges above 
about half the beam charge. These cross sections are 
corrected for interactions in the detectors and interven- 



ing materials, and for secondary interactions in the tar- 
get; they do not require acceptance corrections. Small- 
acceptance charge spectra are used to measure cross sec- 
tions for lighter fragments, and these do require accep- 
tance corrections that monotonically increase with the 
charge change (AZ). 

The small-acceptance charge spectra are found to in- 
variably contain peaks in the region of Z w 3.5. Wo 
have reported these peaks previously; they appear to be 
events in which three helium fragments are detected si- 
multaneously. Here, by comparing cross sections for this 
peak at different beam energies and acceptances, we find 
that the data are best explained by production of 8 Be, 
which instantly decays into two 4 He, in association with 
a third, independently-produced He fragment. 

We presented cross sections for three categories of 
events that are either partially or entirely due to helium 
fragments, which are copiously produced in these reac- 
tions. When the measured cross sections are combined in 
a multiplicity- weighted fashion to yield upper- and lower- 
bound estimates for total helium production cross sec- 
tions, agreement with PHITS and LAQGSM predictions 
is good. The data and both models agree within a factor 
of 2 in almost all cases. Good agreement is seen when 
PHITS is compared to the N and O beam data, and when 
LAQGSM is compared to the Ne and Mg beam data. 

Four hundred eight fragment production cross sec- 
tions presented here were used for comparisons to the 
three models using the Cumulative Relative Uncertainty 
method. All three models yield similar results, as seen 
in Figure 14; the models can all be said to be reasonably 
accurate (U(xi) < 0.5) for 75% of the fragment produc- 
tion cross sections. Overall, PHITS gives the best overall 
agreement, followed by NUCFRG2 and LAQGSM. To 
the limited extent that comparable data are available, 
the present data are in reasonable agreement with earlier 
measurements. 
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